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Abstract. 

We characterize the binary population in the young and nearby OB association Scorpius 0B2 (Sco OB2) using available obser- 
^ ' vations of visual, spectroscopic, and astrometric binaries with intermediate-mass primaries. We take into account observational 

^ : biases by comparing the observations with simulated observations of model associations. The available data indicate a large 

binary fraction (> 70% with 3cr confidence), with a large probability that all intermediate mass stars in Sco 0B2 are part of 
a binary system. The binary systems have a mass ratio distribution of the form fq(q) oc q^i, with y^, x -0.4. Sco OB2 has 
a semi-major axis distribution of the form oc a^" with -/„ a; -1.0 (Opik's law), in the range 5 Ro ^ ^ ^ 5 x 10* Rq. 
' The log-normal period distribution of Duquennoy & Mayor (1991) results in too few spectroscopic binaries, even if the model 

, binary fraction is 100%. Sco 0B2 is a young association with a low stellar density; its current population is expected to be 

very similar to the primordial population. The fact that practically all stars in Sco 0B2 are part of a binary (or multiple) system 
demonstrates that multiplicity is a fundamental factor in the star formation process, at least for intermediate mass stars. 
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1 . Introduction In this paper we aim to recover the properties of the popula- 

tion of binaries that result from the formation process: the pri- 



Over the past decades observations have indicated that a large 
fraction of stars are part of a binary or multiple system. 
Apparently, multiplicity is an important aspect of the star for- 
mation process. Binaries also play a vital role in explaining 
many spectacular phenomena in astrophysics, e.g. supernovae 



mordial binary population, which is defined as the population 
of binaries as established just after the gas has been removed 
from the forming system, i.e., when the stars can no longer ac- 
crete gas from their surroundings (Kouwenhoven et al. 2005). 

^ , „ ,. . ,„„„ ,..„ , , „ XT- The dynamical evolution of stars of a newly born stellar pop- 
type la (Yungelson & Livio 1998; HiUebrandt & Niemeyer w • • a j i. c ^.v. u 
■^^^ •' ulation IS influenced by the presence of gas. Alter the gas has 
2000), short and long gamma-ray bursts (Fryer et al. 1999), , j i r r.i, £ . 

" " ^ V J been removed 1 - 2 Myr after the formation of the nrst massive 

OB runaway stars (Blaauw 1961; Hoogerwerf et al. 2001; , i.- w • i xr . j i. . n i 

■' ' o stars, the binary population is only affected by stellar evolution 

Gualandris et al. 2004), and binary systems with compact rem- , - t- ■ ^ ■ ^ c ■ 

■' ■' ^ and pure N-body dynamics. From a numerical pomt of view, 

nants such as X-ray binaries, millisecond pulsars, and double , j j i. j 

, ^ „ . the primordial binary population can be considered as a bound- 



neutron stars (van den Heuvel 1994; Fryer & Kalogera 1997) 
Binary systems are also known to strongly affect the dynamical 
evolution of dense stellar clusters (Hut et al. 1992; Portegies 
Zwart et al. 2001, 2007). This is an excellent motivation to 
characterize the outcome of the star forming process in terms 



ary between hydrodynamical simulations and N-body simula- 
tions. Hydrodynamical simulations of a contracting gas cloud 
(e.g.. Bate et al. 2003; Martel et al. 2006) produce stars. After 
the gas is removed by accretion and the stellar winds of the 
, . ,. . ,,. most massive O stars, pure N-body simulations (e.g., Portegies 

of multiplicity and bmary parameters. ^wart et al. 2001; Kroupa et al. 1999; Portegies Zwart et al. 

Send offprint requests to: M.B.N. Kouwenhoven e-mail: 2007) can be used to study the subsequent evolution of star 
t.kouwenhovenOsheffield. ac.uk clusters and the binary population. 
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OB associations are well suited for studying the primor- 
dial binary population. They are young, with ages ranging from 
~ 2 Myr (Ori OB lb, Cyg OB2; see Brown et al. 1999; Hanson 
2003) to ~ 50 Myr (a Persei, Cep OB6; see de Zeeuw et al. 
1999). Among OB associations younger than i 20 Myr only 
a handful of the most massive systems have changed due to 
stellar evolution. The effects of dynamical evolution are ex- 
pected to be limited due to their young age and low stellar den- 
sity (< 0.1 Mq pc"^). Moreover, OB associations cover the full 
range of stellar masses (e.g. Brown 2001), in contrast to the 
T associations, their low-mass counterparts. 

In this paper we focus on Scorpius Centaurus (Sco OB2), 
the nearest young OB association, and thus a prime candidate 
for studying the binary population. The proximity of Sco OB2 
(118 - 145 pc) facihtates observations, and the young age 
(5 - 20 Myr) ensures that dynamical evolution has not signif- 
icantly altered the primordial binary population since the mo- 
ment of gas removal. The membership and stellar content of 
the association were established by de Zeeuw et al. (1999) us- 
ing Hipparcos parallaxes and proper motions, and its binary 
population is relatively well- studied. 

Due to selection effects, it is not possible to observe the bi- 
nary population in Sco 0B2 directly. The dataset is hampered 
by the selection of the targets and instrumental limits on the 
observable range in semi-major axis, period, eccentricity and 
mass ratio. The fact that the observed binary population is bi- 
ased makes it difficult to draw conclusions about the true binary 
population. However, by using the method of simulating ob- 
servations of modeled stellar populations (e.g. Kouwenhoven 
2006), it is possible to put constraints on the binary population. 
We accurately model the selection effects of the six major bina- 
rity surveys of Sco 0B2, and compare simulated observations 
with the true observations, to determine the properties of the 
current binary population in Sco OB2. 

The organization of this paper is as follows. In § 2 we 
briefly describe the method and terminology that we use to re- 
cover the true binary population. In § 3 we discuss the Sco OB2 
association and membership issues. In § 4 we describe the 
available datasets with visual, spectroscopic, and astrometric 
binaries in Sco OB2, and outline our models for the respective 
selection effects. In § 5-7 we recover the mass ratio distribution 
and the semi-major axis distribution, and we constrain the ec- 
centricity distribution for the binary population in Sco OB2, 
respectively. Conclusions on the intrinsic binary fraction of 
Sco OB2 are drawn in § 8. The possible differences between 
the current binary population and primordial binary population 
in Sco OB2 are discussed in § 9. In § 10 we compare our re- 
sults with those of others, and we discuss the vahdity of our 
assumptions. Finally, we summarize our main results in § 1 1 . 

2. Method and terminology 

We recover the binary population in Sco OB2 from observa- 
tions using the method of simulating observations of modeled 
stellar populations. This method is extensively described in 
Kouwenhoven (2006), and is briefly summarized below. 

With increasing computer power, it has become possible 
to create sophisticated models of star clusters and OB asso- 



ciations. One can compare these simulated associations with 
the observations of real associations in order to constrain the 
properties of the intrinsic binary population. However, this can- 
not be done directly, as the interpretation of the observational 
dataset is hampered by selection effects. Only a small (and 
biased) subset of the binary population is known. With the 
method of simulating observations of simulations (S.O.S.) one 
characterizes the selection effects, and applies these to the sim- 
ulated association. The simulated observations that are then ob- 
tained can be compared directly with the real observations (see, 
e.g., Kouwenhoven 2006). 

In order to recover the binary population in Sco OB2 
we simulate OB association models with different properties. 
We compare each model with observational data, by simu- 
lating observations for each major binarity survey. With this 
comparison we identify which association model is consistent 
with the observations, and thus constrain the binary popula- 
tion in Sco OB2. Kouwenhoven (2006) shows that this is a safe 
method to derive the binary population. As long as the param- 
eter space (of the binary population) is fully searched, and as 
long as the selection effects are weU-modeled, this method al- 
lows recovery of the intrinsic binary population, as well as the 
uncertainties on each derived property. Furthermore, unlike the 
S.O.S. method used in this paper, the traditional method of cor- 
recting for selection effects (using a "correction factor") may 
lead to erroneous or unphysical results. 

In this paper we make several assumptions when recovering 
the binary population in Sco OB2. In our model we consider 
only single stars and binary systems; no higher order multiples 
are assumed to be present. In § 10. 1 we will briefly return to the 
consequences of this assumption. We assume the distributions 
of the different observed parameters to be independent of each 
other: 

/bp(Mi, M2, a, e, i, m, Q, M) 

= fM„M2(MuM2)faia)feie)Mi)Uoj)faiQ.)fMiM) , 

(1) 

where M\ and M2 are the primary and companion mass, a 
the semi-major axis, e the eccentricity, ( the inclination, o) 
the argument of periastron, Q the position angle of the line 
of nodes, and Ai the mean anomaly at some instant of time. 
Alternatively, one can replace faia) by the orbital period distri- 
bution fpiP). 

In our models the overall binary fraction for the associa- 
tion can be described with a single number, independent of the 
primary mass Fm(Mi) = Fm- Observations of other stellar pop- 
ulations have suggested that the binary fraction may depend on 
primary mass or system mass, and that binary fraction tends 
to increase with increasing primary mass (e.g. Preibisch et al. 
1999; Sterzik & Durisen 2004), though a quantitative descrip- 
tion for FyiiMi) is still unavailable. In our models we therefore 
adopt a binary fraction independent of primary mass, so as to 
keep our description for the binary population in Sco OB2 as 
simple as possible. Note that Kouwenhoven (2006) has shown 
that selection effects may introduce a trend between binary 
fraction and primary mass in the observations, even though 
such an intrinsic trend may not be present in reality. However, 
our assumption does not influence the results significantly, as 
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in this paper we only study the population of binary stars with 
an intermediate mass primary, mostly of spectral type B or A. 
As this corresponds to a small mass range, we neglect the pos- 
sible correlation between binary fraction and primary mass (see 
§ 10.4 for a further discussion). 

For the same reason, we assume the semi-major axis a and 
eccentricity e to be independent of primary mass. The inde- 
pendence of the semi-major axis a with respect to the eccen- 
tricity e may be a good approximation, as observations suggest 
that these parameters are only mildly correlated for solar-type 
stars in the solar neighbourhood (e.g., Duquennoy & Mayor 
1991; Heacox 1997); see also § 7. Note, however, that even 
in the case that this dependence is absent in the intrinsic pop- 
ulation, a correlation may still be present in the observations 
due to selection effects. We assume that the inclination i, the 
argument of periastron co, the position angle of the ascend- 
ing node Q., and the mean anomaly M at some instant of 
time are independent of each other and of all other parame- 
ters. Finally, we assume that the binary systems have a random 
orientation in space (which is not necessarily implied by the 
previous assumptions). Even in the unhkely case that binary 
systems do not have a random orientation, the results do not 
change measurably (see Kouwenhoven 2006). Note that the pri- 
mary and companion mass distributions are never independent, 
/m,,M2(Mi, M2) /m, (Mi)/m2(M2), as by definition Mi > M2. 

For reasons of simplicity, we ignore the interaction between 
close binary stars; our models do not include Roche Lobe over- 
flow or conmion envelope evolution. Low-mass contact bina- 
ries, such as cataclysmic variables, WUMa binaries, and sym- 
biotic stars generally appear on a timescale which is signifi- 
cantly longer than the age of Sco OB2. The higher-mass con- 
tact binaries, such as high-mass X-ray binaries, close binaries 
with mass reversal, and double pulsars could be present, or may 
have escaped the association as runaways. Due to the youth of 
Sco OB2, a few of the closest binaries may have evolved into 
such objects. The non-inclusion of this close binaries, however, 
is unlikely to afl'ect our conclusions on the primordial binary 
population, as we adopt a lower limit to the period of 12 hours 
(§ 6.1). If these binaries are present, our inferred binary fraction 
(§ 8) may be slightly underestimated. 

Throughout this paper we denote the (intrinsic) probability 
density function of a binary parameter x as /^(x) and its cumu- 
lative distribution as F^dx). The corresponding observed distri- 
butions for a binarity survey are denoted as fxix) and F^ix), 
respectively. 



3. The Sco 0B2 association 

Sco OB2 is currently the best studied OB association. It 
consists of three subgroups: Upper Scorpius (US), Upper 
Centaurus Lupus (UCL) and Lower Centaurus Crux (LCC) 
(e.g., Blaauw 1964; de Zeeuw et al. 1999). These three sub- 
groups are likely the result of triggered star formation (e.g., 
Blaauw 1991; Preibisch & Zinnecker 1999, 2007), and in turn 
may have triggered star formation in the p Ophiuchus region. 
Several properties of the three subgroups of Sco OB2 are listed 
in Table 1. 



Preibisch et al. (2002) performed an extensive study of 
the single star population of the US subgroup of Sco 0B2. 
They combine their observations of PMS-stars with those of 
Preibisch & Zinnecker (1999) and de Zeeuw et al. (1999) 
and derive an empirical mass distribution in the mass range 
0.1 Mo < M < 20 Mo (Equation 2). Lodieu et al. (2006) 
on the other hand studied the low-mass and substellar popula- 
tion of Sco OB2 and find a best-fitting value a - -0.6 ±0.1 
of the mass distribution fmiM) oc M" in the mass range 
0.01 - 0.3 Mq. The results of both studies overlap in the re- 
gion 0.1 - 0.3 Mq. In this region Preibisch et al. (2002) find 
a slope a = -0.9 + 0.2 of the mass distribution, while Lodieu 
et al. (2006) find a = -0.6 ±0.1. The slight difference between 
the measured slopes is likely statistical. It is clear, however, that 
the mass distribution for Sco OB2 has a break at a certain value 

in (or near) the mass range 0.1 - 0.3 M©. For this reason we 
model the mass distribution fuiM) in Sco OB2 as follows: 



fM{M) = 



dN(M) 
dM 



^-0.6+0.1 

yy^-0.9±0.2 
j^-2.S±0.5 
^-2.6±0.5 



^min 

0.6 Mo 
2 Mo 

< 



<M <Mp 
<M < 0.6 Mo 

< M < 2 Mo ' 

< M < 20 Mo 

(2) 

<: 0.3 Mo. Note 



where Mmm ^ 0.01 Mq and 0.1 ^ 
that our adopted prescription for fuiM) roughly corresponds 
to the mass distribution derived by Kroupa (2001), while it is 
slightly steeper than Salpeter (a = -2.35) in the intermediate- 
mass regime. 

In our analysis we focus on deriving the properties of the 
intermediate mass binary population, as ample observations of 
these are available. Due to a lack of systematic surveys for 
binarity among low-mass stars in Sco 0B2 we cannot con- 
strain these. For these reasons, the form of the mass distribu- 
tion fhiiM) for M 1.5 Mq is irrelevant, unless both the pri- 
mary and the companion are directly drawn from fM{M). In 
§ 5 we show that the observations exclude the latter possibility, 
given any reasonable value of M^i^ and Mp. In § 5 we will also 
show that binary systems in Sco OB2 are well described with a 
primary mass distribution fmiM) and a mass ratio distribution 
fqiq), so that the exact values of and Mp are irrelevant. 

3.1. The model for Sco 0B2 

We create association models using the STARLAB simulation 
package (see, e.g., Portegies Zwart et al. 2001). The properties 
of the stellar and binary population are projected onto the space 
of observables using an extension of the STARLAB package. 
We adopt a Plummer model (Plummer 1911) with a projected 
half -mass radius of 20 pc, and assume virial equilibrium. Note 
that, as in this paper we do not evolve the models over time, the 
latter assumptions do not afi'ect our results. 

In our model for Sco OB2 we adopt a distance of 130 pc 
(the median distance of the confirmed members of Sco OB2) 
and an age of 5 Myr. Although the subgroups UCL and LCC 
are older than US, the systematic error introduced by our choice 
of the age is small. In our models we slightly overestimate the 
luminosity of stars in the UCL and LCC subgroups, but this af- 
fects only the stars close to the detection limit (see § 3.2), and 
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140^ 
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221' 
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65 


19 


4 
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0.61 


0.52 


LCC 


118^ 


-355 


17-23" 


0.055 


180' 
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57 


9 


1 


0.37 


0.56 


0.44 


Model 


130 


20 


5 


0.00 


9000 


varying 


varying 








varying 


varying 


varying 



Table 1. Properties of the subgroups Upper Scorpius (US), Upper Centaurus Lupus (UCL), and Lower Centaurus Crux (LCC) of Sco 0B2, and 
of our model for Sco OB2. Columns 2^ list for each subgroup its distance, effective radius, and age. Column 5 lists the median interstellar 
extinction towards each subgroup. Column 6 lists the number of confirmed Hipparcos members of each subgroup, and is followed by the ob- 
served number of singles, binaries, triples, and higher-order systems among the confirmed members in columns 7-10, taken from Kouwenhoven 
et al. (2007). Finally, columns 11-13 Ust the observed multipUcity fraction, non-single star fraction, and companion star fraction among the 
confirmed members (see Kouwenhoven et al. 2005, for a definition of these fractions). Note that the latter quantities are lower limits due to the 
presence of unresolved binary and multiple systems. In the bottom row we list the properties of our Sco OB2 model. The number of systems 
N = S + B (i.e., singles and binaries) used in our model includes substellar objects with masses down to 0.02 Mq. References: (1) de Zeeuw 
et al. (1999), (2) de Geus et al. (1989), (3) Preibisch et al. (2002), (4) Mamajek et al. (2002), (5) de Bruijne (1999). 



does not affect the properties of our simulated observations sig- 
nificantly. The error in the age neither affects the interpretation 
of the observed mass ratio distribution, as each observed mass 
and mass ratio is derived from the absolute magnitude of the 
stars, assuming the correct age for the subgroup, and using the 
Hipparcos parallax for each star individually. The distribution 
of these observed mass ratios are then compared with those of 
the model. 

We adopt the extended Preibisch mass distribution in 
Equation 2 for our model of Sco OB2. We make the assump- 
tions that (1) the mass distribution for the subgroups UCL 
and LCC is identical to that of US, (2) we adopt a minimum 
mass Mmin = 0.02 Mq (i.e., we do not consider planetary- 
mass objects in our mass distribution), and (3) we adopt = 
0.1 Mq for the mass distribution. As in our study we focus on 
intermediate-mass binaries, assumptions (2) and (3) are only 
of importance if both components are directly drawn from the 
mass distribution. In § 5 we discuss this issue and we wiU 
show that the exact values of Mmin and are irrelevant for 
our study. 

Preibisch et al. (2002) estimate that the US subgroup con- 
tains approximately 2525 single/primary stars in the mass 
range 0.1 - 20 M©. With the extension to lower mass in 
Equation 2 the number of singles/primaries is higher, as we 
also include the very low mass stars and brown dwarfs. For a 
minimum mass M^in = 0.01 M© and a value Mp - OA Mq, and 
assuming that the UCL and LCC subgroups have an equal num- 
ber of singles/primaries, the total number of singles/primaries 
in Sco OB2 is approximately 9 200. For a value =s 0.1 Mq 
in Equation 2 the number of singles/primaries is approximately 
8 000. Free-floating planets (Mi ^ 0.02 Mq) are not included in 
the above statistics. We will therefore adopt N = S +B = 9 000 
systems (Mi > 0.08 Mq) in our simulations, where S is the 
number of single stars, and B the number of binary systems. 

3.2. Photometry 

We obtain the magnitude of each simulated star in the opti- 
cal and near-infrared bands using the isochrones described in 
Kouwenhoven et al. (2005). These isochrones consist of mod- 
els from Chabrier et al. (2000) for 0.02 Mq < M < 1 Mq, 



Palla & Stabler (1999) for 1 Mo < M < 2 Mq, and Girardi 
et al. (2002) for M > 2 Mq. We adopt the isochrone corre- 
sponding to an age of 5 Myr and solar metallicity. By adopting 
5 Myr isochrones we overestimate the brightness of 20 Myr 
old stars by ~ 0.05 mag in JHKs for stars with M ^ 1 M© 
and by ~ 0.5 mag in JHKs for stars with M < 1 Mq. The 
error introduced by the metalhcity (~ 0.05 mag in JHKs) is 
negligible for our purposes: see Kouwenhoven et al. (2007) for 
a more detailed description of these matters. The Hipparcos 
magnitude Hp for each star is derived from its V magnitude and 
V - I colour, using the tabulated values listed in the Hipparcos 
Catalogue (ESA 1997, Vol. 1, § 14.2). For each star we con- 
vert the absolute magnitude into the apparent magnitude using 
the Hipparcos parallax of each star. We do not include inter- 
stellar extinction in our models. Sco OB2 is practicaUy cleared 
of gas. The median visual extinction for the member stars of 
the three subgroups is Avxss = 0.5 mag, Avxscl = 0.06 mag, 
and, AvLCC = 0.05 mag, respectively (de Bruijne 1999), which 
translate to values of J us ~ 0.05 mag, A^Ts.ucl ~ 0.006 mag, 
and, AksXCC ~ 0.006 mag in the near-infrared (Mathis 1990). 
For the purpose of our study the interstellar extinction can thus 
be neglected, in particular for the study of the near-infrared 
surveys of Shatsky & Tokovinin (2002), Kouwenhoven et al. 
(2005), and Kouwenhoven et al. (2007). 

3.3. Sco 0B2 membership 

De Zeeuw et al. (1999) have published a census of the stellar 
content and membership of nearby (^ 1 kpc) OB associations. 
They present a list of 521 members of the Sco OB2 associa- 
tion, based on the Hipparcos position, proper motion, and par- 
allax of each star. Of these members, 120 are in the US sub- 
group, 221 in UCL, and 180 in LCC. Due to the Hipparcos 
completeness limit, most of the confirmed members are bright 
(y ^ 8 mag) and mostly of spectral type B, A, and F. 

In the analysis of the observational data (§ 4) we consider 
only the confirmed members of Sco OB2 (i.e., those identified 
by de Zeeuw et al. 1999), all of which are in the Hipparcos 
catalogue. Among the stars observed by Hipparcos it is un- 
hkely that a Sco OB2 member star is not identified as such. 
On the other hand, it is possible that non-members are falsely 
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the proportion P as a function of V with three tine segments: 
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Fig. 1. The completeness of the Hipparcos catalogue in the Sco OB2 
region, as a function of V magnitude. The diamonds represent the ra- 
tio between the number of stars in the Hipparcos catalogue and the 
number of stars in the TYCHO-2 catalogue, in each V magnitude bin. 
The comparison above is made for the Sco 0B2 region, and is similar 
for each of the three subgroups of Sco OB2. The solid line represents 
the model for the completeness adopted in this paper (Equation 3). 



classified as members of Sco OB2; the so-called interlopers. 
The fraction of interlopers among the "confirmed" Sco OB2 
members stars is estimated to be ~ 6% for B stars, ~ 13% for 
A stars, and ~ 22% for F and G stars (see Tables A2 and CI in 
de Zeeuw et al. 1999). The interlopers among B and A stars are 
Ukely Gould Belt stars, which have a distance and age compa- 
rable to that of the nearby OB associations. In our analysis we 
assume that all confirmed members in the list of de Zeeuw et al. 
(1999) are truly member stars, and do not attempt to correct for 
the presence of interlopers. 

The Hipparcos completeness limit is studied in detail by 
Soderhjehn (2000). His prescription for the completeness is 
based on all entries in the Hipparcos catalogue. However, many 
OB associations were studied in detail by Hipparcos, based 
on candidate membership lists. Due to the Hipparcos crowd- 
ing limit of 3 stars per square degree, only a selected subset 
of the candidate members of Sco OB2 was observed (see de 
Zeeuw et al. 1999, for details), which significantly complicates 
the modeling of the Hipparcos completeness. We therefore cal- 
ibrate the completeness of Hipparcos in the Sco OB2 region 
by comparing the number of Hipparcos entries with the num- 
ber of stars of a given magnitude in the same region. We use 
the TYCHO-2 catalogue for this comparison. The TYCHO-2 
catalogue is complete to much fainter stars than Hipparcos. 
In Figure 1 we show the proportion P of stars that is in the 
Hipparcos catalogue, relative to the number of stars in the 
TYCHO-2 catalogue, as a function of V magnitude. We model 



I for y< 6.80 mag 

logP(y) = I 1.18 - 0.17 V for 6.80 mag < V < 8.88 mag . 
[ 9.18 - 1.07 V for 8.88 mag < V 

(3) 

The TYCHO-2 catalogue is 99% complete down to V = 11 mag 
and 90% complete down to V = 1 1.5. The completeness of the 
Hipparcos catalogue for V ^ 11 mag is therefore not accu- 
rately described by Equation 3. However, this does not affect 
our results, as the surveys under study only include the bright- 
est members of Sco OB2. In this model we ignore the fact that 
the Hipparcos completeness also depends on spectral type. 

Apart from the large membership study of de Zeeuw et al. 
(1999), several others have been performed. Several recent 
studies have focused on the membership of low-mass objects, 
in particular of brown dwarfs (e.g. Kraus et al. 2005; Slesnick 
et al. 2006; Lodieu et al. 2006). These studies often focus on 
a small group of suspected members in a specific region of 
Sco OB2, and the membership is mostly based on the (less ac- 
curate) photometric method. As no census on membership of 
low-mass stars and brown dwarfs is currently available, we do 
not include these in our analysis. 

4. Observations of binary systems in Sco 0B2 

A large fraction of the Sco OB2 member stars is known to be 
part of a binary or multiple system. In Table 1 we have pro- 
vided an overview of the observed binary fraction in the as- 
sociation, for which we included all known binary and mul- 
tiple systems in Sco OB2. In total there are 266 known com- 
panions among the 521 confirmed members of Sco OB2, most 
of which are intermediate mass stars. The multiplicity fraction 
in Sco OB2 is at least 40% among these stars, assuming that 
all proposed companions are indeed physical companions. The 
references for these binary and multiple systems are listed in 
Table 2. Figure 2 shows several observed parameter distribu- 
tions. These observed distributions are not representative of the 
intrinsic distributions, as selection effects prohibit the detection 
of a significant fraction of the companion stars. Furthermore, 
it is possible that several of the reported companions are spu- 
rious, including for example bright background stars that are 
projected close to a Sco OB2 member star. 

The known binary systems in Sco OB2 were discovered by 
different observers, using various techniques and instruments. 
As each of these observing runs is characterized by specific 
selection effects, it is difficult to study each of these in detail. 
We therefore focus primarily on a subset of the surveys: those 
of Kouwenhoven et al. (2005, 2007); Shatsky & Tokovinin 
(2002); Levato et al. (1987); Brown & Verschueren (1997), and 
those detected by Hipparcos (ESA 1997). We refer to these pa- 
pers and the corresponding datasets as K05, K06, SHT, LEV, 
BRV, and HIP, hereafter. An overview of the number of ob- 
served targets and detected binary systems in each dataset is 
presented in Table 3. Combined, these datasets contain a large 
fraction of the known binary and multiple systems in Sco OB2. 
The selection effects for each of these datasets can be modeled, 
making it possible to use the method of simulated observations. 
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Fig. 2. Properties of the observed binary population in Sco 0B2. Only the 521 confirmed members of Sco 0B2 are considered. The top-left 
panel shows the angular separation distribution for visual binaries, at the moment of observation. The other panels show the distribution over 
radial velocity amplitude Ki, the projected semi-major axis Ui sin/, the eccentricity e, the mass function f(M), and the period P, derived for 
the orbits of the specttoscopic binaries (SBl and SB2), and for HIP78918, the only astrometric binary in Sco OB2 with an orbital solution. The 
measurements shown in this figure include those of multiple systems. Above each panel we indicate the number of companions for which the 
corresponding orbital element is available. Spectroscopic and astrometric binaries without an orbital solution are not included. 



Reference 


Detection method 


rtiCIlCdl CL al. \Z,\J\jj ) 


opccuoacupic 


rtiiuciocii CI al. \i.yyj) 




od-icgti CI til. 1 yy-r j 


V l&Utll 


Barbier-Brossat et al. (1994) 


Spectroscopic 


Batten etal. (1997) 


Spectroscopic 


Buscombe & Kennedy (1962) 


Spectroscopic 


Chen et al. (2006) 


Visual 


Couteau (1995) 


Combination 


The Double Star Library 


Combination 


Duflotetal. (1995) 


Specttoscopic 


Hartkopf et al. (2001) 


Visual 


Jilinski et al. (2006) 


Specttoscopic 


Jordietal. (1997) 


Eclipsing 


The Hipparcos and TYCHO Catalogues 


Astromettic 


Kraicheva et al. (1989) 


Specttoscopic 


Kouwenhoven et al. (2005) 


Visual 


Kouwenhoven et al. (2007) 


Visual 


Lindroos (1985) 


Visual 


Malkov (1993) 


Combination 


Mason (1995) 


Visual 


McAhster et al. (1993) 


Visual 


Miscellaneous, e.g. SIMBAD 


Combination 


Miuraetal. (1992) 


Visual 


Nitschelm (2004) 


Specttoscopic 


Pedoussaut et al. (1996) 


Specttoscopic 


Shatsky & Tokovinin (2002) 


Visual 


Sowell& Wilson (1993) 


Visual 


Svechnikov & Bessonova (1984) 


Combination 


Tokovinin (1997) 


Combination 


Worley & Douglass (1997) 


Combination 



Table 2. References to literature data with specttoscopic, asttomet- 
ric, eclipsing, and visual binaries among the Hipparcos members of 
Sco 0B2. The data for a number of binary systems in Sco 0B2 is taken 
from several catalogues. This table is similar to the one presented in 
Kouwenhoven et al. (2005), but is updated with recent discoveries. 



In the following sections we describe these five datasets, and 
discuss our model for the selection effects. 

In our approach we follow the modeling of the selection 
effects for visual, spectroscopic, and astrometric binarity sur- 
veys, which is extensively discussed in Kouwenhoven (2006). 



A summary of the modeled selection effects for each of the sur- 
veys is given in Table 4. With the sample bias we refer to the 
process of selecting the targets of interest. A selected sample 
for a survey usually consists of a group of stars with specific 
properties, such as the solar-type stars in the solar neighbour- 
hood, or the B stars in an association. In the case of a binary 
survey among the members of an OB association, the observer 
may erroneously include a bright background star, assuming 
that it is an association member. We consider this part of the 
sample bias, although we do not discuss this aspect in this pa- 
per. When the selected targets are surveyed for binarity, obser- 
vational constraints are responsible for the instrument bias. For 
example, the minimum and maximum detectable angular sep- 
aration of binary stars is determined by the properties of the 
telescope and the detector. We include in the instrument bias 
the selection effects imposed by the telescope-instrument com- 
bination and atmospheric conditions. We additionally include 
the bias that results from the difficulties of identifying compan- 
ions. For example, faint companions at a large separation of 
their primary may not always be identified as such, due to the 
confusion with background stars. 



4.1. K05 — Kouwenhoven etal. (2005) observations 

Kouwenhoven et al. (2005) performed a near-infrared adap- 
tive optics binarity survey among A and late-B members 
of Sco OB2. Their observations were obtained with the 
ADONIS/SHARPII-i- system on the ESO 3.6 meter telescope 
at La Silla, Chile. Adaptive optics was used to obtain high spa- 
tial resolution, in order to bridge the gap between the known 
close spectroscopic and wide visual binaries. The survey was 
performed in the near-infrared, as in this wavelength regime 
the contrast between the components of a binary system with 
a high mass ratio is less than in the visual regime. All tar- 
gets were observed in the Ks band, and several additionally 
in the J and H bands. K05 selected their sample of A and 
late-B targets from the list of confirmed Hipparcos members 
that were identified by de Zeeuw et al. (1999). All targets have 
6 mag £V £9 mag, which corresponds to similar limits in the 
Ks band. 
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Abbreviation 


Reference 


Dataset 




^orig 


Abused 


^used 


K05 


Kouwenhoven et al. (2005) 


Visual 


199 


74 


199 


60 


K06 


Kouwenhoven et al. (2007) 


Visual 


22 


29 


22 


18 


SHT 


Shatsky & Tokovinin (2002) 


Visual 


115 


25 


87 


23 


LEV 


Levato et al. (1987) 


Spectroscopic 


81 


61 


53 


39 


BRV 


Brown & Verschueren (1997) 


Spectroscopic 


156 


91 


71 


47 


fflP 


Ferryman et al. (1997) 


AsttometricA'isual 


521 


125 


521 


125 



Table 3. An overview of the datasets used to derive the properties of the binary population in Sco OB2. Columns 1-3 list the dataset acronym, 
the reference, and the type of binary studied in the dataset. Columns 4 and 5 list the number of targets in the original dataset, and the number 
of companions found for these targets. Columns 6 and 7 list the number of targets and companions used in our analysis. This dataset is smaller 
than the original dataset, as we do not include the non-members of Sco OB2 in our analysis and at most one companion per targeted star in 
the case of a multiple systems. The datasets partially overlap, which is taken into account when these are combined in the following sections. 
We list in this table the total number of spectroscopic binaries, including the radial velocity variables (RVVs; irrespective of their true nature), 
SB Is, and SB2s. For the Hipparcos observations we list the number of entries in the categories (X), (O), (G), (C), and (S), among the confirmed 
members of Sco OB2. 



With their observations K05 are sensitive to companions 
as faint as Ks ~ 15.5 mag, corresponding to the brightness of 
a massive planet in Sco OB2. Due to the large probabihty of 
finding faint background stars in the field of view, K05 clas- 
sify all secondaries with Ks > 12 mag as background stars, 
and those with Ks < 12 mag as candidate companion stars. The 
Ks = 12 criterion separates companion stars and background 
stars in a statistical manner, and is based on the background 
star study of SHT. A member of Sco OB2 with Ks - 12 mag 
has a mass close to the hydrogen-burning Umit. The follow- 
up study of K06 with VLT/NACO (see § 4.2) has shown that 
the Ks - 12 criterion correctly classifies secondaries as com- 
panions in 80 - 85% of the cases. With their survey K05 
find 151 secondaries around the 199 target stars. Out of these 
151 secondaries, 74 are candidate companions {Ks < 12 mag), 
and 77 are background stars {Ks > 12 mag). K05 find that 
the mass ratio distribution fq{q) for late-B and A type stars in 
Sco OB2 is consistent with fq{q) oc q'""-^^, and exclude random 
pairing between primary and companion. 

AAA. Treatment of the K05 dataset 

All 199 targets in the K05 dataset are confirmed members 
of Sco OB2, and are therefore included in our analysis. We 
use in our analysis a subset of the companions identified in 
K05. Several targets in the ADONIS survey have more than 
one candidate or confirmed companion. In this paper we do 
not study triples and higher-order multiples; we consider at 
most one companion per target star. For each of these candi- 
date multiple systems we include the (candidate) companion 
that is most likely a physical companion. For HIP52357 we in- 
clude the companion with {p,Ks) = (0.53", 7.65 mag), as it is 
brighter and closer to the target star than the candidate compan- 
ion with {p,Ks) - (10.04", 11.45 mag). For the same reason, 
we do not include the wide and faint candidate companion of 
HIP61796 with (p, Kg) = (12.38", 1 1.86 mag) in our analysis. 
K05 find two bright and close companions of HIP76001, with 
{p,Ks) = (0.25", 7.80 mag) and {p,Ks) = (1.48", 8.20 mag), 
respectively. Although HIP76001 is likely a physical triple, 
we choose to retain only the innermost candidate companion. 
K05 find a bright secondary separated 1.8" from HIP63204. 



With their follow-up study, K06 find a close companion at 
p - 0.15". K06 show that this close companion is physical, 
while the secondary at 1.8" is optical; we do not consider the 
latter secondary in our analysis. 

HIP68532 and HIP69113 are both confirmed triple sys- 
tems, each with a primary and a "double companion". For both 
HIP68532 and HIP691 13, the two stars in the "double compan- 
ion" have a similar separation and position angle with respect 
to the primary, and a similar magnitude. In physical terms, the 
double companions of HIP68532 and HIP69113 could have 
originated from a more massive companion that fragmented 
into a binary. We therefore model the double companions of 
these stars as single companions, taking the average p and (f, 
the combined Ks magnitude, and the total mass of each double 
companion. 

For the comparison with the simulated observations the tar- 
gets HIP77315 and HIP77317 are both considered as individ- 
ual, single stars. The star HIP77317 is known to be a compan- 
ion of HIP77315 at p = 37.37", and is for that reason listed 
as such in K05. This binary system is far too wide to be de- 
tected with the observing strategy of K05; both stars are there- 
fore treated as individual stars. With the ADONIS survey K05 
find three candidate companions of HIPS 1972. Of these three, 
only the secondary at separation 5.04" is a confirmed compan- 
ion in the follow-up study of Kouwenhoven et al. (2007). As 
HIP81972 is near the Galactic plane, the other two secondaries 
are Ukely background stars, and are thus not included in the 
dataset. 

K05 separated the secondaries into candidate companions 
and background stars using the Ks magnitude of each sec- 
ondary. The follow-up study of K06, using multi-colour anal- 
ysis, has shown that several of these candidate companions are 
background stars. We do not consider in our analysis these 
secondaries, indicated with HIP53701-1 {Ks = 8.9 mag), 
HIP60851-1 {Ks = 11.5 mag), HIP60851-2 {Ks = 11.3 mag), 
HIP80142-1 {Ks = 9.51 mag), and HIP80474-1 {Ks = 10.8 
mag) in K06. 

The resulting K05 dataset that we use for our analysis con- 
tains data for 199 targets with a total of 60 companion stars. For 
each of these targets and their companions we use the measure- 
ments given in K05, unless more recent (and more accurate) 
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K05 — (Kouwenhoven et al. 2005) — Visual binaries 



Observer's choice 

Brightness constraint 
Separation constraint 
Contrast constraint 
Confusion constraint 



A and late-B members of Sco OB2 

(incl. Hipparcos completeness) 

5.3 mag < Vi < 9.5 mag, Mi > 1.4 Mq 

Equation 7 

Equation 4 

Ks2 < 12 mag 



Adaptive optics, -50% detection limit 
I I 1 1 1 1 1 1' I I I I 1 1 I I I I 1 1 1 1 1 1 I I I I I I I I I I 1 1 1 1 



K06 — (Kouwenhoven et al. 2007) — Visual binaries 



Observer's choice 

Brightness constraint 
Separation constraint 
Contrast constraint 
Confusion constraint 



A selection (1 1%) of the K05 sample 
5.3 mag <Vi< 9.5 mag. Mi > 1.4 Mq 
Equation 10 
Equation 8 
Not applicable 



SHT — Shatsky & Tokovinin (2002) — Visual binaries 



Observer's choice 

Brightness constraint 
Separation constraint 

(idem) 
Contrast constraint 

(idem) 
Confusion constraint 



B members of Sco OB2 
(incl. Hipparcos completeness) 
Vi < 7.0 mag. Mi > 3.5 Mo 
Non-coronographic mode: Equation 11 
Coronographic mode: Equation 14 
Non-coronographic mode: Equation 13 
Coronographic mode: Equation 15 
A's_2 < 12 mag, J2 < 13 mag, 
and J 2 - Ks2 < 1.7 mag 



LEV — Levato et al. (1987) — Spectroscopic binaries 



Observer's choice 

Brightness constraint 

Contrast constraint 
Amplitude constraint 
Temporal constraint 
Aliasing constraint 
Sampling consttaint 



B members of Sco 0B2 

(incl. Hipparcos completeness) 

Vcomb < 8.1 mag. Ml >3Mo 

Not applicable 

Spectroscopic bias model SB-W 
with T = 2.74 year, Ar = 0.38 year, 
and eRv = 3.1 kms~'. 

Not applied 



BRV — Brown & Verschueren (1997) — Spectroscopic binaries 



Observer's choice 

Brightness constraint 
Contrast constraint 
Amplitude constraint 
Temporal constraint 
Aliasing constraint 
Sampling constraint 



B members of Sco OB2 

(incl. Hipparcos completeness) 

Vi < 7.0 mag, Mi > 3.5 Mq 

Not applicable 

Spectroscopic bias model SB-W, 
with T = 2.25 year, Ar = 0.75 year, 
and eRv = 1.4 km s"'. 

Not applied 



HIP — Hipparcos mission — Astrometric binaries 
Brightness constraint Hipparcos completeness 



AmpUtude constraint 
Temporal constraint 

Aliasing constraint 
Sampling constraint 



Classification into the 
categories (X), (C), (O), or (G) 
depending on the observables 
of each binary system (see Table 6) 



Table 4. An overview of the models for the selection effects used to 
generate simulated observations of simulated OB associations, for the 
six major datasets discussed in Sections 4.1 to 4.6. The sample bias, 
resulting from the choice of the sample alone, includes the observer's 
choice and the brightness constraint. All other constraints result from 
the properties of the telescope, detector, atmospheric conditions, and 
coirfusion with background stars, and are in this paper referred to as 
the instrument bias. For a detailed description of the constraints men- 
tioned in this table we refer to § 4.5 of Kouwenhoven (2006). 



measurements for these stars are presented in the follow-up 
study of K06. For each of the targets HIP63204, HIP73937, 
and HIP79771 a new close companion is resolved by K06, 




12 3 4 

Angular separation (arcsec) 

Fig. 3. The 50% detection limit AKg as a function of angular sep- 
aration p, for the K05 observations (short-dashed curve), the K06 
observations (dotted curve), the non-coronographic SHT observa- 
tions (dash-dotted curve), the coronographic SHT observations (long- 
dashed curve), and the combined SHT observations (solid curve). For 
the K05 and K06 observations the curves represent those for average 
Sttehl ratios of 30% and 24%, respectively. 



which was unresolved in the observations of K05. For these 
three targets we use the properties of the primary star as pro- 
vided by K06. 

The mass of each primary and companion is derived from 
the near-infrared magnitude. If available, the mass of each star 
is taken from K06, who use the JHKs magnitude to derive 
the mass. In all other cases the mass is taken from K05, who 
derive the mass from the Ks magnitude only. The more re- 
cently determined masses allow us to better constrain the prop- 
erties of the binary population in Sco OB2. Finally, the mass 
ratio q = M2/M1 is calculated for each binary system. In the 
Appendix (Table A. 1) we list the properties of the binaries used 
for comparison with simulated observations. 

4.1 .2. Modeling the observational bias of K05 

We model the sample bias in K05 as follows. The authors 
selected the A and late-B members of Sco OB2. As these 
members were identified as such in the Hipparcos membership 
study of de Zeeuw et al. (1999) we first impose the Hipparcos 
completeness (see § 3.3) on the simulated association. Based 
on the properties of the target list of K05, we model the ob- 
server's choice and brightness constraint by removing all tar- 
gets (i.e. singles and primaries) with V < 5.3 mag, all targets 
with V > 9.5 mag, and all targets with M < 1.4 M© from the 
sample. 

We model the detection limit of the K05 observations using 
the analysis presented in K06, who describe these in detail. We 
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study the 50% detection limit (in terms of the magnitude differ- 
ence AKs = K2- Ki) and find its dependence on angular sep- 
aration and Strehl ratio (SR). We parametrize the dependence 
of the detection Umit A^rs,det(p) in magnitudes on Strehl ratio 
(SR) as 



A/srs,det(p) = 

ro 

(22.0-3.75 s(SR))(jD-0.1") 
8.8 + s(SR)(p-2") 



P <Plim,A 
Plim,A < P < 0.5" 

0.5" < p < 2" 
2" <p 



(4) 

where piim.A = 0.2" is the angular resolution of the K05 obser- 
vations. Following the properties of the K05 observations, we 
model i(SR) with 



s(SR) = 2.54- 3.85 X SR. 



(5) 



As an example we plot the detection limit A/rs,det(p) in Figure 3 
for observations with SR = 30%. We simulate the distribution 
over SR by drawing for each target the SR randomly from the 
observed distribution ^r(SR), which is approximated with 



/sr(SR) oc exp 



(SR-z^sr ) 



2^ 



5% < SR < 50% , (6) 



where jUsr = 30% and ctsr = 5%. 

K05 considered only the secondaries with Ks <\2 mag as 
physical companions. The foUow-up study of K06 has shown 
that this Ks = 12 criterion indeed correctly classifies most of 
the companions and background stars. We therefore identify 
in our simulated observations only the companions with Ks < 
12 mag as true companions. 

Each measurement is assigned a detection probability 
Z)a(p) as a function of angular separation p. This detection 
probability refers solely to whether or not a companion is pro- 
jected into the field of view. As the field of view is not circular, 
the detection probabihty of a companion is a function of an- 
gular separation. For ADONIS we have a square field of view 
sized 12.76" x 12.76". K05 observed each target four times, 
each time with the target in another quadrant of the field of 
view, so that the efl'ective field of view is La = 5 • 12.76" = 
19.14". The probability that a secondary with an angular 
separation p is in the field-of-view is then given by: 

1 for p < La/2 

Da{p) = ] 1 - (4/7r) arccos(LA/2p) for La/2 < p < La/ V2 , 
for LkI^<p 

(7) 

which is visualized in Figure 4. 

4.2. KOS — Kouwenhoven et al. (2007) observations 

The results of the ADONIS binarity survey performed by K05 

raised several questions, in particular on the absence of faint 
secondaries in the 1 " - 4" separation range, and on the valid- 
ity of the Ks = 12 criterion that K05 used to separate sec- 
ondaries into companion stars and background stars. Although 
SHT and KOS argue that the latter criterion statistically clas- 
sifies the background stars correctly, the correct classification 
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Angular separation (arcsec) 

Fig. 4. In the three imaging surveys for binarity discussed in this pa- 
per (KOS, K06, and SHT) the field-of-view is non-circular. Whether a 
secondary is in the field of view, depends therefore not only on its sep- 
aration p, but also on its position angle ip. This figure shows the prob- 
ability that a secondary is in the field-of-view, as a function of p, as- 
suming random orientation of the binary systems, for the KOS obser- 
vations (short-dashed curve), the K06 observations (dotted curve), the 
non-coronographic SHT observations (dash-dotted curve), the corono- 
graphic SHT observations (long-dashed curve), and the combined 
SHT observations (solid curve). Whether a secondary is detected or 
not, depends additionally on its brightness and on the brightness dif- 
ference with the primary (see Figure 3). 



of the companion stars with Ks ~ 12 mag was still uncertain. 
To address this issue, K06 performed follow-up multi-colour 
JHKs observations of a subset of the ADONIS targets. With 
multi-colour observations, each secondary can be placed in the 
colour-magnitude diagram, and compared with the isochrone 
of the Sco OB2 subgroups. Companion stars are expected to 
be near the isochrone, while background stars are (generally) 
expected to be far from the isochrone. 

The observations described in K06 were carried out 
with the adaptive optics instrument NAOS-CONICA (NACO), 
mounted on the ESO Very Large Telescope on Paranal, Chile. 
A subset of 22 (out of 199) K05 targets were selected for 
follow-up observations. The subset was not randomly selected, 
but preference was given to faint and close background stars, 
to secondaries with Ks ~ 12 mag, and to newly discovered 
candidate companions. K06 analyzed the JHKs observations 
of these 22 stars observed with NACO, including the multi- 
colour ADONIS observations of 9 targets. With their observa- 
tions K06 found three new close companions (of mP63204, 
H1P73937, and H1P79771) that were unresolved in the survey 
of K05. 
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4.2.1 . Treatment of the K06 dataset 

For our analysis we consider all 22 targets observed by K06, all 
of which are confirmed members of Sco OB2. The 9 ADONIS 
targets that were also studied in K06 are not considered here, 
simply because they were not observed in the campaign of 
K06. Around the 22 NACO targets K06 find 62 secondaries, 
of which they classify 18 as confirmed companions (c), 11 as 
possible companions (?), and 33 as background stars (b). 

In our analysis we use the data for 15 (out of 18) con- 
firmed companions, and 5 (out of 8) candidate companions. 
Both HIP68532 and HIP69113 have a tight "double compan- 
ion". We treat each of these as a single companion, by combin- 
ing the separation and mass of the individual companions (see 
§ 4.1.1). The targets fflP67260, HIP79771, and fflP81949 all 
have three (candidate) companions, for which we only include 
the inner companion in our analysis. We do not include the very 
faint secondary of HIP80142, as this is likely a background 
star. For HIP81972 we only include the companion HIP81972- 
3, which is by far the most massive companion, in our analysis. 

The final K06 dataset used in this paper consists of 22 tar- 
gets with 18 companions. Note that when this dataset is 
compared with simulated observations, a discrepancy may be 
present, as the sample was composed to study candidate com- 
panions and background stars with particular properties. In 
Table A. 2 we list the properties of the binaries used for com- 
parison with simulated observations. 

4.2.2. Modeling the observational bias of K06 

A subset of 22/199 = 11% of the targets in the K05 sample 
are observed with NACO by K06. We model the K06 sam- 
ple by randomly drawing 11% of the targets in the simulated 
K05 target sample. Note that in reaUty, the subset was not ran- 
dom (see above); instead, the targets were selected based on 
the properties of their secondaries. The simulated K06 obser- 
vations therefore cannot be directly compared with the results 
of the K06 observations. However, they can be used to find the 
expected number of close and/or faint companions with K06; 
companions that could not be found with the K05 survey. 

We use the 50% detection limit from the analysis presented 
in K06, and parameterize it with the Strehl ratio (SR) of the 
observations. The 50% detection limit as a function of p, for 
targets with a difl'erent brightness is derived using simulations 
(Kouwenhoven et al. 2007, § 3). From the observational data 
we derive a detection limit A/rs,det(p) in magnitude, consisting 
of four Une segments: 

^Ks4ct(p) = 

P < Plim,N 

(2.5 B(SR) - 2.63) (p - 0.1") pu^.N < p < 0.5" 
£(SR)-h0.70(p-2") 0.5"<p<2" 
B(SR) p > 2" 

(8) 

where piijn,N = 0. 1" is the angular resolution of the K06 obser- 
vations and 



B(SR) = 6.86 -I- 8.37 x SR 



(9) 



is the magnitude difl'erence of the faintest detectable source, for 
a given Strehl ratio SR. As an example we show the detection 
limit A/rs,det(p) in Figure 3 for observations with SR = 24%. 
We simulate the distribution over SR by drawing for each tar- 
get the SR randomly from the observed distribution /sr(SR), 
which we approximate with Equation 6, with fisR = 24% and 
csR = 7%. The field of view for the observations of K06 is 
14" X 14". As the field of view is non-circular, the detection 
limit is a function of both angular separation p and position 
angle (p. For our simulated observations, each measurement is 
assigned a detection probability Dn(p) as a function of angular 
separation p, given by 

■ 1 p < Ln/2 

£)n(p) = < 1 - (4/7r) arccos(LN/2p) Ln/2 < p < L^/ V2 , 
( Ln/ < p 

(10) 

where Ln = 14" is the linear size of the field-of-view. 

4.3. SHT — Shatsky& Tokovinin (2002) observations 

Shatsky & Tokovinin (2002) performed an imaging survey for 
binarity among 115 B type stars in the Sco OB2 region. Their 
observations were carried out in 2000 with the near-infrared 
adaptive optics instrument ADONIS at the ESO 3.6 meter tele- 
scope on La Silla, Chile. Their sample is based on the study 
of Brown & Verschueren (1997); see § 4.5. Among the 115 
B-type stars surveyed by SHT, 87 are confirmed members of 
Sco OB2 according to de Zeeuw et al. (1999). Among the total 
sample of these 115 stars SHT find 96 secondaries in the angu- 
lar separation range 0.3" - 6.4", of which they identify 10 as 
new physical companions. The authors conclude that the mass 
ratio distribution fqiq) for B-type stars in Sco OB2 is consistent 
with fq{q) oc q~^-^, and that random pairing can be excluded. 

4.3.1 . Treatment of the SHT dataset 

Near the 87 confirmed members of Sco OB2 targeted by 
SHT, 80 secondaries are found, of which 61 likely optical and 
19 likely physical companions. Of this set of 19 physical com- 
panions, we use a subset of 17 for our analysis. The target 
HD 132200 is probably a physical triple system. As we con- 
sider in this paper only single and binary systems, we do not 
include the widest and faintest component of HD 132200, and 
retain the component with p - 0. 128" and Ks - 5.46 mag. The 
secondary HD133937P is incorrectly reported in SHT. For this 
secondary p = 0.57" and ./ - Ks = 2.06 mag (N. Shatsky & 
A. Tokovinin, private communication). Due to its large J - Ks 
value, and as 7 > 13 mag, this secondary is likely a background 
star. We therefore do not consider HD133937P in our analysis. 

Several targets were not included in the observed sam- 
ple of SHT. These targets were known to have close compan- 
ions and thus not suitable for wavefront sensing. These unob- 
served targets were included in the analysis of SHT though. 
Seven of these (HIP53701, HIP57851, HIP62322, HIP64425, 
HIP741 17, HIP76371, and HIP77840) are confirmed members 
of Sco OB2. We include these non-observed targets in our 
analysis, either as single or as a binary system, depending on 
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whether their companions would have been detected with the 
SHT observing strategy. Technically, the non-inclusion of a set 
of stars falls under the "observer's choice". The latter con- 
straint is difficult to model, as it would involve modeUng of 
pre-SHT observations of close binaries, as well as the determi- 
nation whether or not such a binary is suitable for wavefront 
sensing. We choose, however, to manually add these stars to 
the list of observed targets, as the properties of these stars and 
their companions are well-understood (making detailed models 
of the observer's choice redundant). Note that non-inclusion of 
these unobserved stars introduces a small bias, as these stars 
would have been surveyed by SHT if their companions were 
unknown at that time. The member stars fflP5785 1 , HIP62322, 
HIP741 17, HIP76371, and H1P77840 were reported as visually 
resolved (C)-binaries in the Hipparcos catalogue. We use the 
angular separation and magnitude of these components as given 
in the catalogue, and include the stars in the sample. HIP64425 
is a known triple system (Tokovinin 1997) for which we use the 
massive inner binary in our analysis. We treat the non-observed 
star HIP53701 as a single star, as K06 have shown that its sec- 
ondary is a background star. 

For the stars observed by SHT, we derive the mass of tar- 
get and companion star from the Ks magnitude, using the evo- 
lutionary models described in Section 3. For the stars that 
are analysed by SHT, but not observed by these authors (see 
above), we derive the mass using the V band magnitude and 
Hipparcos Hp magnitude. For each star we adopt the distance 
given by the Hipparcos parallax and the age of the subgroup of 
which the target is a member (Table 1). 

The final dataset from the SHT survey that we use in our 
analysis of the binary population in Sco OB2 comprises 87 tar- 
gets with 23 physical companions. The properties of these 
23 companions are Usted in Table A.3. 

4.3.2. Modeling the observational bias of SHT 

We require that all targets are confirmed Hipparcos members 
of Sco OB2, and so first impose the Hipparcos detection limit 
on the simulated association. We model the brightness con- 
straint of the SHT observations by adopting a minimum mass 
of 3.5 Mo, and a minimum brightness of V = 7 mag for the 
targets. 

SHT show the typical detection limit of their observations 
in their Figure 3. The detection limit is obviously different for 
the observations with and without the coronograph. The ob- 
servations with coronograph are deeper, and the observations 
without the coronograph provide a larger range in angular sep- 
aration. A companion star is detected if it is observed either in 
the coronographic mode or in the non-coronographic mode. 

In the non-coronographic observations, each companion is 
assigned a detection probabiUty D^cip) as a function of its 
separation p. In the non-coronographic observations, SHT ob- 
served each target twice in the non-coronographic mode, both 
times with the target in a quadrant of the detector. Due to the 
square shape of the detector, and due to the observing strategy, 
the position angle is of importance to whether a companion at 
separation p is in the field-of-view. We model this dependence 



by assigning a probability Dnc(p) that a companion is in the 
field of view, depending on p. For a square field-of-view of 
a detector with hnear size L, and a separation K between the 
two observations (along the diagonal of the field-of-view), the 
probabiUty is given by 



1 



+ - arcsm 



+ - arcsm 



=■ arccos 



(?) 
(?) 



p<ri 

ri < r2l ^/2 

r2^<p<r^ 

r3 <p<r2 

r2<p 



(11) 



where 



r, = \Lyl2-\K, r2=\Lyl2 + \K, r^ = ^1(^2 + ^2). 

(12) 

For the non-coronographic observations of SHT, the linear di- 
mension of the detector is L = 12.76", and the translation along 
the diagonal of the field of view is K = 8.5" (see Figure 1 in 
SHT, for details). We model the detection limit of the observa- 
tions without the coronograph with 

(0 forp<piim,s 
A/s:s,det = { 8.32 logp + 6.83 forpu^.s <p< 1.46" , 
[8.2 forp>1.46" 

(13) 

where pum.s =0.1" is the angular resolution of the SHT obser- 
vations (based on Figure 3 in SHT). 

SHT additionally observe each target using the corono- 
graph. They do not perform their coronographic observations 
in mosaic-mode; only one pointing is used. Each measurement 
is therefore assigned a detection probability Dcip) as a function 
of angular separation p, given by 



Dc(p) 



for p < dc 

1 for < p < L/2 

1 - ^ arccos (^) for L/2 < p < L/ V2 
for L/y/2<p 



(14) 



where L - 12.76" and dc - 1" is the radius of the corono- 
graph. Based on Figure 3 in SHT, we model the detection limit 
of the observations with the coronograph with 



A/s:s,det = ] 8.321ogp-H6.83 
10.3 



for p < dc 

for dc<p< 2.62" . (15) 
for p > 2.62 



Finally, we combine the simulated observations in corono- 
graphic and non-coronographic mode. We consider a binary 
system as detected, if it is observed in at least one of the two 
modes. 

SHT additionally studied the background star population in 
the Sco OB2 region. Due to the large number of background 
stars, it is likely that a very faint or red secondary is a back- 
ground star. SHT classify a secondary as a background star if 
7 > 13 mag, if ^5 > 12 mag, or if J - Ks > 1.7 mag (unless 
the secondary is a known companion). In our model for the se- 
lection eff'ects, we adopt these limits in magnitude and colour 
when obtaining the simulated observations. 
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4.4. LEV — Levato etal. (1987) observations 

Levato et al. (1987) performed a large radial velocity sur- 
vey for binarity among early-type stars in the Sco OB2 re- 
gion. They performed their observations in May 1974 with the 
0.9 meter and 1.5 meter CTIO telescopes, and in 1976 with 
the 2.1 meter telescope at KPNO. Their sample consists of 

81 candidate members of Sco OB2, and is based on that of 
Slettebak (1968) who composed a list of suspected Sco OB2 
members for a study on stellar rotation. All except 4 of the 

82 targets of Slettebak (1968), and 3 additional targets were 
observed by LEV. The spectral type of the observed targets 
ranges from BOV to AOV. The targets in the sample have 
2.5 mag < V < 8.1 mag. 

On average, each star is observed over an interval of (T) = 
2.74 year, with a spread of ctt = 0.68 year. Each target is ob- 
served 5-12 times, with an average observing interval {AT} = 
0.38 year and a corresponding spread of cr^r = 0.14 year. 
For each target LEV list the internal error in the radial ve- 
locity measurements. Averaged over all targets, this error is 
(fRv) = 3.1 kms \ with aspreadof 1.0 kms"'; approximately 
90% of the targets have crv > 2 km s"^ . 

In their Table 3 LEV list their conclusions on binarity. Of 
the 53 confirmed members of Sco OB2 that they observed, 
14 have a constant radial velocity (within the measurement er- 
rors), 23 have a variable radial velocity (RVV), 8 are SB 1, and 
8 are SB2. Given these observations, the spectroscopic binary 
fraction is at least (8 -i- 8)/53 = 30%, in the case that all RVVs 
are spurious. If aU reported RVV targets are indeed binaries, the 
observed spectroscopic binary fraction is (8-i-8+23)/53 = 74%. 
The observed spectroscopic binary fraction is a lower limit for 
Fu, as binaries that are unresolved in the survey of LEV (e.g. 
visual binaries) are not included in these statistics. 

4.4.1 . Treatment of the LEV dataset 

For the comparison between the observational data and the sim- 
ulated observations, we only consider the 53 confirmed mem- 
bers (according to de Zeeuw et al. 1999) of Sco OB2 that 
LEV observed. In their Tables 3 and 4, LEV include the star 
HIP76945 (HD 140008), a confirmed member of the UCL sub- 
group. LEV did not observe this SB2, but take the orbital ele- 
ments from Thackeray & Hutchings (1965). Our simulations 
indicate that the radial velocity variability of a binary with 
properties such as those of HIP76945 would practically always 
be detected in a survey similar to that of LEV. It is unclear, 
however, whether LEV would have been able to derive the or- 
bital elements for this binary, i.e., if they would have detected 
it as an SBl or SB2. We include HIP76945 as an SB2 in the 
dataset, as Thackeray & Hutchings (1965) were able to derive 
the orbital elements several decades before the study of LEV. 
For a subset of the targets LEV derive the orbital elements. In 
their Table 4, LEV list the elements of 22 targets, of which 
16 ais confirmed members of Sco OB2. In Table A.4 we list 
the properties of these 16 SBl and SB2 systems from the LEV 
dataset that are confirmed members of Sco OB2. We also list 
the 23 radial velocity variables (RVVs), for which the orbital 



elements are unavailable. The LEV dataset consists of 53 tar- 
gets, of which 16 -I- 23 = 39 are detected as binary systems. 

4.4.2. Modeling the observational bias of LEV 

In this paper we consider only the confirmed Hipparcos mem- 
bers of Sco OB2, i.e., first impose the Hipparcos detection limit 
on the association. We model the choice of the sample of LEV 
by removing all binary systems with a combined magnitude 
fainter than V = 8.1 mag from the simulated observations. 

We model the instrument bias of LEV using windowed 
sampling (SB-W; see Kouwenhoven 2006). Briefly summa- 
rized; in order to compare the model predictions with the ob- 
servations, we simulate the detection of the spectroscopic bina- 
ries in our models as follows. We obtain radial velocity mea- 
surements of all binary systems in the simulated association, at 
regular intervals AT" for a time-span T (windowed sampling). 
We assume a value for the measurement error e^y, which is 
constant over the time of observations. If the radial velocity 
measurements show a spread significantly larger than the error, 
the binarity is detected. For each single star and binary system 
we test the hypothesis that the observed velocity measurements 
{v,) result from a constant velocity. We calculate thsx^ of the 
set of A'v radial velocity measurements: 

where v is the mean of the measurements {v,}. We then calculate 
the probability p that;^'^ is drawn from the;k'^-distribution: 

/, = l-r(iy,i;r2), (17) 

where v = A^v - 1 is the number of degrees of freedom. High 
values (p « 1) indicate that our hypothesis (that the radial ve- 
locity is constant) is true, and that the measurements are likely 
the result of statistical noise. Values of p close to zero indicate 
that the observed variations in the radial velocity are real. We 
classify objects with radial velocity sets with p < 0.0027 (cor- 
responding to the 3 cr confidence level) as binary systems, while 
the other targets are marked as single stars. In our modeUng 
of spectroscopic binaries we thus only determine whether a bi- 
nary is detected or not; we do not discriminate between spectro- 
scopic binaries of type SBl, SB2 or RVV. ModeUng the latter 
difference is sophisticated and depends on a significant num- 
ber of parameters. For example, the Nyquist theorem requires 
that at least two measurements should be obtained per orbital 
period. Additional constraints are imposed by the properties of 
the binary system: the spectral type of the star (in particular the 
number of spectral Unes), the brightness of the system, and the 
values of the radial velocity amplitude Ki (relative to eRv), the 
eccentricity e and the argument of periastron oj. It is not trivial 
to model these. 

In our model for the LEV observations we use windowed 
sampling, adopting an observing run of T = 2.74 year, an ob- 
serving interval AT = 0.38 year, and a radial velocity accuracy 
of = 3.1 kms"'. The latter assumption is a simplification, 
as the value of erv is slightly different for each observation in 
the LEV dataset (with a spread of ~ 1 kms~^). A star is more 
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easily detected if e^y < 3.1 km s"\ and less easy if is larger. 
Our simulations show, however, that our assumption of a con- 
stant Crv introduces an error significantly smaller than the error 
introduced by low-number statistics, justifying our assumption. 

4.5. BRV— Brown & Verschueren (1997) 
observations 

Brown & Verschueren (1997) studied stellar rotation among 
members of the Sco OB2 association. The observations were 
carried out between 1991 and 1993 using the ECHELEC spec- 
trograph at the ESO 1.52 meter telescope on La Silla, Chile. 
The sample of BRV contains the pic-Hipparcos candidate and 
established members of Sco OB2, based on the studies of 
Blaauw (1964), Bertiau (1958), and de Geus et al. (1989). 
The observations and data reduction procedure are described 
in detail in Verschueren et al. (1997), and the results on du- 
plicity are described in Verschueren et al. (1996) and Brown 
& Verschueren (1997). Their sample consists of 156 targets in 
the Sco OB2 region, mostly of spectral type B. They find that 
~ 60% of the binary systems exhibit a significant radial ve- 
locity variation. After combination of their data with those of 
LEV and those of the Bright Star Catalogue (Hoffleit & Jaschek 
1982; Hoffleit et al. 1983), they obtain a binary fraction of 74%. 

4.5.1 . Treatment of the BRV dataset 

Among the 156 observed targets there are 71 confirmed mem- 
bers of Sco OB2 (1 8 in US, 30 in UCL, and 23 in LCC). Among 
these 71 targets, 7 are SBl, 10 are SB2, 30 are RVV, and 
12 have a constant radial velocity (CON). For 12 targets, insuf- 
ficient measurements are available to make a statement about 
the radial velocity variation. 

Two out of the 30 RVV binaries are known to exhibit radial 
velocity variation due to line profile variability. HD120324 is a 
non-radial pulsator and HD 136298 is a yS Cephei variable. For 
both stars, this is likely the reason that they are classified as 
RVV. These stars are therefore not considered as binary sys- 
tems in our analysis. 

The final BRV dataset used in our analysis consists of 
71 confirmed members. Of these targets, 12 are spectroscop- 
ically single, 7 are SBl, 10 are SB2, 28 are RVV, and 12 have 
insufficient data to determine whether the radial velocity is 
variable. The binary fraction is thus at least (7 + 10)/71 ^ 
24%, if none of the RVV and CON targets are binary, and 
(71 - 12)/71 = 83% if all RVV and CON targets are bi- 
nary. Among the target stars with sufficient data to make a 
statement on duplicity (i.e., wiffiout the 12 CON targets in- 
cluded), ffie observed binary fraction is 17/59 = 28% at least 
and 45/59 = 76% at most. 

4.5.2. Modeling the observational bias of BRV 

We model the choice of ffie BRV sample in a way similar as 
we did for the LEV dataset. Each target is observed three times 
over an interval of T = 2.25 year, so ffiat AT - 0.75 year. 
Following ffie reduction of the original data (Verschueren et al. 



1996), we classify each target wiffi a radial velocity variation 
larger than 3eRv = 4.2 km s"> as a RVV (see § 4.4.2). 

4.6. HIP — Hipparcos observations 

In the Hipparcos catalogue, the (possible) binary nature of a 
target is indicated in field H59 wiffi ffie flags (X), (O), (G), (C), 
or (V), and in field H61 wiffi ffie flag (S). For ffie targets wiffi a 
(C) flag both stars in the binary system are resolved, but no or- 
bital motion is detected. These systems are considered as visual 
(or optical) binaries. Hipparcos entries with an (X)-flag have a 
stochastic solution. These stars exhibit an apparent motion sig- 
nificantly larger ffian the statistical uncertainties, although no 
double star solution could be found. For entries wiffi an (O)- 
flag, at least one of the orbital elements could be derived from 
the apparent motion. Entries with a (G)-flag show a significant 
acceleration in ffie apparent motion, but no solution could be 
found. These are likely long-period binaries. The (V)-flagged 
entries are variability-induced movers. For this group of bina- 
ries, ffie photocenter exhibits apparent motion due to variability 
of one of the components. Finally, (S)-flagged entries are sus- 
pected non-single stars. These targets are effectively single as 
observed by Hipparcos (ESA 1997, § 2.1), alffiough no con- 
vincing single-star astrometric model could be fit to the obser- 
vations. Several of ffie (S) binaries are also (X)-flagged, indi- 
cating ffiat these are likely non-single. 

Among the 521 confirmed members of Sco OB2, 46 are 
candidate or confirmed astrometric binaries — i.e., ffiose in 
ffie categories (X), (O), (G), and (S) — in the Hipparcos cat- 
alogue: 12 in US, 17 in UCL, and 17 in LCC. An additional 
79 Sco OB2 members are classified as (C)-binaries; these are 
visuaUy resolved binaries. Table 5 lists for ffie three subgroups 
of Sco OB2 the number of entries in each of the Hipparcos cat- 
egories. As for ffie binaries in the category (C) no orbital mo- 
tion is detected, we will consider this group as visual binaries. 
(V)-binaries are not present among the confirmed Sco OB2 
members. HIP78918 is the only member of Sco OB2 wiffi 
an orbital solution (O). The three (X)-flagged members of 
Sco OB2 are also flagged as suspected non-single stars (S). 

4.6.1 . Treatment of the Hipparcos dataset 

In our analysis we consider each target in the categories (X), 
(O), (G), and (C) in Table 5 as a binary system. Binarity among 
the stars in the (S) category (the "suspected non-single" tar- 
gets) is raffier uncertain. We therefore compare our results wiffi 
ffie Hipparcos data, with and wiffiout the suspected (S) binaries 
included. Note that not all targets in the category (X) are nec- 
essarily binary systems. For example, the flag (X) of HIP80763 
(or Sco) may be due to ffie extended nature of the star, which is 
surrounded by a dust-shell (Cruzalebes et al. 1998). This may 
induce photocentric motion ffiat is not related to binarity. 

4.6.2. Modeling the observational bias of Hipparcos 

The observer's choice and sample bias for the Hipparcos 
member list of Sco OB2 are discussed in Section 3.3; our 
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Subgroup 




(X) 


(O) 


(G) 


(C) 


(S) 


^^M,X0G 


^MJCOGS 




US 


120 


1 





4 


15 


8(1) 


4.2 ± 1.5% 


10.0 ± 2.1% 


12.5 ± 3.0% 


UCL 


221 





1 


9 


36 


7(0) 


4.5 ±1.5% 


7.7 ± 1.9% 


16.3 ± 3.0% 


LCC 


180 


2 





6 


28 


11(2) 


4.4 ± 1.6% 


9.4 ± 2.4% 


15.6 ± 3.0% 


Sco OB2 


521 


3 


1 


19 


79 


26(3) 


4.4 ± 1.0% 


8.8 + 1.4% 


15.2 ± 1.8% 



Table 5. Candidate and confirmed astrometric binaries in the Hipparcos catalogue. For each subgroup we list the number Af* of known members, 
the number of stochastic (X), orbital (O), acceleration (G), component (C), and suspected (S) binaries in the Hipparcos catalogue. For each 
(S) binary we list between brackets how many of these are also (X)-flagged. The last three columns list the "astrometric binary fraction" — 
including the (X), (O), (G) binaries — without the (S) binaries and with the (S) binaries included, and the Hipparcos "visual" binary fraction, 
for the (C) binaries only. (V) binaries are not present in Sco OB2. 



Constraints on p and AH,, 


Period constraints 


Solution 


Symbol 


Elements 


2 < <p) < 100 mas or AH^ > 4 


P < 0.1 year 


Stochastic 


(X) 


no 


2 < <p) < 100 mas or AHp > 4 


0.1 <P< 10 year 


Orbital 


(O) 


yes 


2 < p < 100 mas or AHp > 4 


5 < P < 30 year 


Acceleration 


(G) 


no 


0.1 < p < 100 arcsec and AHp < 4 


P > 30 year 


Resolved 


(C) 


no 


Not modeled 




Suspected 


(S) 


no 


Not modeled 




VIM 


(V) 


no 



Table 6. A model for the instrument bias of the Hipparcos catalogue, based on the analysis of Lindegren et al. (1997). The binary systems 
satisfying the above constraints are resolved with Hipparcos in our models. For the comparison between the observations and the simulated 
observations, we consider two sets of Hipparcos binaries: the visual binaries and the astrometric binaries. No orbital motion is detected for the 
(C) binaries; these are visually resolved and therefore technically visual binaries. The Hipparcos astrometric binaries contain the targets with 
(X), (O), (G), and optionally (S) entries. No difference between the latter categories is made for the comparison with the astrometric binaries. 
Binary systems that do not satisfy the constraints listed in this table remain undetected in our simulated observations for Hipparcos. We do not 
model the (V)-binaries (variability-induced movers; VlMs) and (S)-binaries (suspected non-single stars). Note that in our model we overpredict 
the number of binaries in categories (X), (O), and (S), as not all binaries with the properties above are detected by Hipparcos as such. 



adopted model to describe these biases is given in Equation 3. 
The instrument bias for Hipparcos was described in detail in 
Lindegren et al. (1997) and is summarized in Table 6 (see 
Kouwenhoven 2006, for further details). 

The binaries in category (C) are considered as visual bina- 
ries, as both components of such binaries are resolved, while 
no astrometric motion is detected. The detection of the (C)- 
binaries is modeled using the prescription in Table 6 and is then 
compared with the observations. 

The binaries in the categories (X), (O), and (G) are con- 
sidered as astrometric binaries, and are modeled using the pre- 
scriptions in Table 6. As we do not model the Hipparcos obser- 
vations in detail, we are unable to accurately predict in which of 
these categories each Hipparcos target falls. The simple model 
that we adopt for the Hipparcos biases results in an overlap 
between the properties of the stars in these categories, and fur- 
thermore, we overpredict the number of stars observed in these 
categories. Lindegren et al. (1997) have analyzed the proper- 
ties of the binary systems in each of the categories. But this 
does not mean that each binary system with these properties 
is observed as such. In our model we make the latter assump- 
tion, resulting in an overestimation of the number of binaries 
detected by Hipparcos, and an overlap between the modeled 
categories (O) and (G). We therefore combine the number of 
detected (X), (O), and (G) binaries in the simulated observa- 
tions, and use the resulting astrometric binary fraction Fm.xog 
as an upper limit for the comparison with the observations. 

We do not include stellar variability in our model, and are 
therefore unable to model the (V)-flagged binaries (field H59). 
We do not model the (S)-flagged binaries (field H61) either, as 



the determination whether a target is classified as such in the 
Hipparcos catalogue is based on an internal and external com- 
parison performed by the two Hipparcos data reduction teams. 

5. Recovering the pairing function and mass ratio 
distribution 

In this and the following sections we discuss the determination 
of the binary population in Sco OB2 from observations. We 
adopt the strategy described in Kouwenhoven (2006). First we 
derive the pairing function and mass ratio distribution in § 5. 
We recover the semi-major axis, the period distribution and the 
binary fraction in § 6, and the eccentricity distribution in § 7. 
Finally, we present a discussion on the derived binary fraction 
and its associated error in § 8. 

Kouwenhoven (2006) discusses five possible ways of pair- 
ing the components of a binary system. These include random 
pairing (RP), primary-constrained random pairing (PCRP), and 
three variants of primary-constrained pairing (PCP-I, PCP-II, 
and PCP-III). For the models with random pairing the primary 
mass Ml and companion mass M2 are both drawn from the 
mass distribution. For PCRP both masses are drawn from the 
mass distribution, with the additional constraint that the com- 
panion mass is smaller than that of the primary. In the three 
PCP pairing models the primary mass is drawn from the mass 
distribution, and the companion mass is derived from the mass 
ratio q = M2/M1 which is drawn from a distribution fqiq). The 
difference between the three PCP models lies in the treatment 
of the very low mass companions that are generated. Each of 
these five pairing functions result in a different binary popula- 
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Fig. 5. How well can the intermediate-mass binary population in Sco 0B2 be described by a mass ratio distribution of the form /^{q) oc q'l'] 
Left: the goodness-of-fit for the comparison between the observed mass ratio distribution (from the K05 and SHT datasets) and that predicted by 
models with different values of -y^. A large value for the goodness-of-fit means that the model predictions are consistent with the observations. 
The Icr, 2a and 3(T confidence limits for rejection of the model are indicated with the horizontal dotted lines. Middle: the goodness-of-fit 
for the comparison between the observed visual binary fraction (from the K05 and SHT datasets) and that predicted by the models with an 
intrinsic binary fraction of 85% (bottom curve), 90%, 95% and 100% (top curve). Right: the predicted (SBl, SB2 and RVV) spectroscopic 
binary fraction (for the LEV dataset) as a function of -y^. The five curves indicate the results for a model binary fraction of 60% (bottom curve), 
70%, 80%, 90% and 100% (top curve). The bottom horizontal dotted line indicates the minimum observed spectroscopic binary fraction, where 
only the SB Is and SB2s are included, assuming that the radial velocity variations in all RVVs are caused by line-profile variability rather 
than binarity. The top horizontal line indicates the maximum spectroscopic binary fraction, assuming that all RVVs are indeed binaries. Taken 
together, the three panels indicate that the mass ratio distribution for intermediate mass stars in Sco OB2 can be described by fq{q) oc q^i with 
exponent ya ~ -0.4 (vertical dotted line in each panel). An intrinsic binary fraction close to 100% is required in order to produce the large 
visual and spectroscopic binary fractions that are observed. See § 5 for a further discussion of this figure. 



tion. Kouwenhoven (2006) has shown that for binary systems 
with intermediate mass primaries the pairing functions PCP-I, 
PCP-II, PCP-III practically give the same results. As we focus 
on intermediate mass binaries in this paper, we will therefore 
consider three pairing functions in our analysis: RP, PCRP, and 
PCP (where PCP represents either PCP-I, PCP-II, or PCP-III). 

As shown by Kouwenhoven (2006) the mass ratio distri- 
bution fqiq) resulting from pairing function RP or PCRP de- 
pends strongly on the mass distribution fM(M), in particular on 
the shape of the mass distribution in the brown dwarf regime. 
The mass distribution for Sco OB2 is fairly well constrained 
(Equation 2), although the exact values of the parameters Mmin 
and Mp are as yet unknown. However, the mass distribution 
can still be used to rule out pairing functions RP and PCRP for 
Sco OB2, by considering the extreme values for Mmin and Myj. 

Shatsky & Tokovinin (2002) and Kouwenhoven et al. 
(2005) have already shown that their observed mass ratio dis- 
tribution is inconsistent with random pairing (RP). Below we 
show that both pairing functions RP and PCRP can be excluded 
based on their observations. The two free parameters in the 
mass distribution of Sco 0B2 (Equation 2) are in the range 
Mrain ^ 0.01 Mo and 0.1 Mq ^ ^ 0.3 Mq. The shape 
of the resulting mass ratio distribution for these pairing func- 
tions depends on the exact values of M^in and M^. Below we 
will show that even for the "most favourable" values of M^in 
and Mj3 both RP and PCRP can be excluded. 

The mass distribution in Equation 2 results in a large num- 
ber of low-mass stars. The probability to obtain a binary con- 
sisting of two intermediate or high-mass stars is small (e.g. 
Kouwenhoven 2006). We use this property of pairing functions 



RP and PCRP to show that these are inconsistent with the ob- 
servations. Let Q denote the ratio between the number of bina- 
ries with mass ratio q > 0.8 and the number of targets A'taigets: 

# binaries with ^ > 0.8 
Q= 7j • (18) 

^targets 

The value of Q increases with increasing M,„in and increasing 
Mj3 due the smaller probability of drawing a low-mass object 
from the mass distribution. Furthermore, Q is proportional to 
the intrinsic binary fraction Fm- For pairing functions RP and 
PCRP, the largest value of Q is therefore reached if Mmin ~ 
0.01 Mo, Mp « 0.3 Mo, and Fm = 100%. 

We simulate two associations with the latter properties (one 
with pairing function RP and one with PCRP) and extract 
the sample of K05 using the constraints listed in Table 4. 
The resulting intrinsic values of Q for the K05 sample are 
Qrp ~ Gpcrp ~ 0.004%. These values are upper limits because 
of the adopted values for Mmin, Mp, and Fm. Also, the corre- 
sponding observed values Q of the simulated observations are 
significantly smaller due to the non-detection of very close and 
wide binaries, binaries with a high mass ratio, etc. 

Among the 199 targets in the K05 sample, 10 binary sys- 
tems with q > 0.8 are detected. The observed fraction of bi- 
naries with q > 0.8 is therefore Qkos = 5.0 ± 1.6%. Due to 
selection effects K05 have certainly missed several binaries 
with q > 0.8 at separations smaller than the spatial resolution, 
or outside the field of view, so that the true value is 2ko5 > 
(2ko5 = 5.0 + 1.6%. As Qrp <s; Qkos and Qpcrp «: Qkos, 
both pairing functions RP and PCRP can thus be excluded with 
strong confidence. 
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The observations are consistent with pairing function PCP, 
which is characterized by a mass ratio distribution fq{q). The 
unknown slope a and the value Mp for the mass distribution 
(Equation 2) are now irrelevant, as the distribution over com- 
panion masses among intermediate mass stars is independent of 
a. The mass ratio distribution and binary fraction for high-mass 
and intermediate mass targets is the same for the three pairing 
functions (Kouwenhoven 2006). As we do not have detailed 
information on binarity among low-mass stars in Sco OB2, we 
cannot discriminate between PCP-I, PCP-II, and PCP-III. A de- 
tailed membership study for low-mass stars in Sco OB2, fol- 
lowed by a detailed binary study, is necessary to establish the 
difference. Kouwenhoven (2006) also shows that for the three 
PCP pairing functions the mass ratio distribution for binaries 
with a high-mass primary is approximately equal to the gener- 
ating mass ratio distribution fq{q) of an association. 

The mass ratio distribution among intermediate mass stars 
in Sco OB2 was obtained by both K05 and SHT SHT find a 
mass ratio distribution of the form fq{q) oc q^i with jq - -0.5 
in their B star survey, and K05 find = -0.33 in their A and 
late-B star survey. Subsequently, we adopt a mass ratio dis- 
tribution of the form fq{q) oc q^i and study for which value 
of jq the simulated observations correspond best to the obser- 
vations of K05 and SHT. We compare the observed mass ra- 
tio distribution fq{q) with that of the simulated observations 
using the Kolmogorov-Smirnov (KS) test, and we compare 
the observed binary fraction Fm with the predictions using 
the Pearson test. In both cases, we test the hypothesis that 
the observations and simulated observations are realizations 
of the same underlying association model. All models have 
a semi-major axis distribution of the form /«(«) oc a"' with 
0.5 R© < a < 5 X lO*" Ro (see § 6.1, § 6.2, and § 6.5 for a 
motivation) and a thermal eccentricity distribution /e(e) = 2e 
with < e < 1 (see § 7 for a motivation). We perform our sim- 
ulations with different intrinsic binary fractions, ranging from 
Fm = 60% to 100%. 

Figure 5 shows the results of this comparison for models 
with a varying value of jq and a varying intrinsic binary frac- 
tion Fm- The left-hand panel shows the probability associated 
with the KS comparison between the observed fq{q) and that 
of the simulated observations. A small value for the goodness- 
of-fit means that the model can be excluded with high confi- 
dence. The three horizontal dotted hues indicate the Icr, 2cr and 
3cr confidence limits for rejection of the model, respectively. 
Models outsize the range -0.6 jq ^ -0.3 can be rejected 
with Icr confidence. The best-fitting models have jq ^ -0.4 
(vertical dotted line), which is the value we adopt for Sco OB2. 
However, jq is also constrained using the additional informa- 
tion provided by the observed visual and spectroscopic binary 
fraction. 

The middle panel of Figure 5 indicates how well the ob- 
served visual binary fraction for Sco OB2 compares with the 
predictions of the models with different values of jq and Fm- 
Models with a small value of jq underpredict Fm. Most bina- 
ries have a small mass ratio in this case, making it more difficult 
to detect the companion star, and hence resulting in a lower vi- 
sual binary fraction. Models with a Fm < 95% can be excluded 



with 2cr confidence or more, for any value of < 0.2 (under 
the condition that our assumptions hold). 

The right-hand panel of Figure 5 shows the predicted spec- 
troscopic binary fraction as a function of jq and Fm. Note 
that, unlike the middle panel of Figure 5, we do not show the 
goodness-of-fit between the observed and predicted spectro- 
scopic binary fraction, as we only have lower and upper lim- 
its of the former (see below). More spectroscopic binaries are 
detected for models with a large jq, as the radial velocity am- 
plitude increases with higher companion mass. The horizontal 
dotted lines in this panel indicate the limits for the observed 
spectroscopic binary fraction: 30% ^ Fm ^ 74% (with Ukely 
Fm ~ 74%). For models with k -0.6, the models with 
an intrinsic binary fraction Fm ~ 100% are most consistent 
with the observations. Note that the spectroscopic binary frac- 
tion as predicted by the models is rather small as compared to 
the observed spectroscopic binary fraction, even for an intrin- 
sic binary fraction Fm = 100%. This apparent underprediction 
is also present in the following sections, where we study dif- 
ferent semi-major axis and period distribution, and may result 
from the presence of triple and higher-order systems among 
the members of Sco OB2. In Section 10.1 we will return to this 
issue. 

A good model for Sco OB2 should correctly predict fqiq), 
the visual binary fraction and spectroscopic binary fraction. 
The distribution fp(p) suggests that -0.6 ^ jq ^ -0.3, with 
a best-fit for jq ^ -0.4, while the observed binary fractions are 
best fitted by models with a large value of y^. A combination 
of these data is not trivial, as it involves a choice for the relative 
weights given to each of the goodness-of-fits for each quantity. 
For example, if fpip) and the visual binary fraction are given 
equal weight, the combined best-fitting value is jq ^ -0.35, 
while if fp(p) is a ten times heavier weight, the combined best- 
fitting value is jq « -0.41. Inclusion of the spectroscopic bi- 
nary fraction is not trivial, as we only have lower and upper 
limits for this quantity. An inspection of the three panels sug- 
gets that the mass ratio distribution can be well described by 
the expression fq(q) oc q'>'i with jq x -0.4. This value is brack- 
eted by, and consistent with the values derived by K05 and 
SHT. The comparison further indicates that the intrinsic binary 
fraction is larger than « 95%, with Icr confidence. 

6. Recovering the semi-major axis distribution, the 
period distribution and the binary fraction 

The orbital size distribution of a binary population can be 
quantified in several ways. Observers of visual binaries of- 
ten express their results with a semi-major axis distribution 
faia), derived from the observed angular separation distribu- 
tion. Observers of spectroscopic binary systems on the other 
hand, often express their results with an orbital period distri- 
bution fp{P), as the period of a spectroscopic binary can often 
be measured directly. Authors of theoretical/simulation papers 
mostly describe the orbital size distribution using the orbital 
energies (e.g. in units of kT). Over the past several decades, ob- 
servational studies have brought forward two widely accepted 
distributions: a flat distribution in logarithmic semi-major axis, 
and a log-normal period distribution. In the this section we 
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compare these distributions with the observations. We wish to 
stress that the discussion of /„(«) and fp(P) should nowhere be 
taken to imply that a power-law distribution in a would result 
in a log-normal distribution in P. 

The flat distribution in log a, commonly known as Opik's 
law, has been derived for a wide range of stellar populations 
(e.g., Opik 1924; van Albada 1968; Vereshchagin et al. 1988; 
Poveda & AUen 2004), and is equivalent to 

fa{a) OC d^' Uain <a< «max , (19) 

with = - 1 . For a set of binaries of total mass Mt, this results 
in an orbital period distribution 

/p(P)ocA4^«^i^/3p(2r.-i)/3^ (20) 

which, in the case of 7a = -1, is a flat distribution in log P (see 
Section 4.D in Kouwenhoven 2006). Throughout this paper we 
will consider the distribution /a(a) and study for which values 
of Ya the models are consistent with the observations. 

Duquennoy & Mayor (1991) studied binarity among solar- 
type stars in the solar neighbourhood and find a log-normal 
period distribution: 

/DM(/^)'^exp|-^i^i^5^| Prmn<P<Pm^, (21) 

here /ip = logP = 4.8, <Tp = cr^ogp = 2.3, and P is in 
days. The latter distribution is often used as the standard ref- 
erence for the orbital size distribution of a binary population. 
This log-normal period distribution results in an approximately 
log-normal semi-major axis disli ibulion, the shape of which is 
mildly dependent on the distribution over binary mass Mj. For 
a set of binaries of total mass Mj, the resulting semi-major axis 
distribution is exactly described by a log-normal distribution 
with mean and width 

loga = f logP - i log (^) and crioga = |o-iogP (22) 

respectively, where a is in astronomical units (see Section 4.D 
in Kouwenhoven 2006). Throughout this paper we will con- 
sider the general form for the log-normal period distribution 
fp;ti,a-{P), and study for which combination of jU and cr the mod- 
els are consistent with the observations. 

In our binarity dataset for Sco OB2, most information on 
the orbital size distribution is provided by the observed angular 
separation distribution fpip), the observed visual binary frac- 
tion and the observed spectroscopic binary fraction. The latter 
two parameters are additionally linearly dependent on the in- 
trinsic binary fraction Fm of Sco OB2, so we derive Fm simul- 
taneously. 

For the comparison of the observed angular separation dis- 
tribution fp(p) we use the angular separation measurements of 
the combined observations of the K05 and SHT datasets (see 
§4.1 and § 4.3). The observed (visual) binary fraction Fyi is 
30 + 4% for the K05 dataset and 26 + 4% for the SHT dataset. 
The visual binary fraction for the combined K05/SHT dataset 
is Fm = 31 ± 4% (taking into account the 23 targets that ap- 
pear in both datasets, and their companions). We will use the 
latter value for comparison with the simulated observations. 
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Table 8. The widest known binary systems in Sco 0B2. For each of 
these binary systems we list the angular separation, the Hipparcos par- 
allax, an estimate for the semi-major axis Cest = ZJtanp, the subgroup, 
and the spectral type of the primary. The last column lists the refer- 
ence. Note that this list must be incomplete, as very wide binaries are 
difficult to detect. Furthermore, several of these binaries may be opti- 
cal due to confusion with background stars, i.e., not physically bound. 
References: (1) Lindroos (1985), (2) Worley (1978), (3) Worley & 
Douglass (1997), (4) Oblak (1978), (5) Tokovinin (1997). 



For the comparison between the observed spectroscopic bi- 
nary fraction and the predictions we use the observations of 
LEV. The LEV sample (see § 4.4) consists of 53 confirmed 
members of Sco OB2. Among these targets there are 8 SB Is, 
8 SB2s, 23 RVVs, and 14 targets with a constant radial ve- 
locity (within the measurement errors). For the 23 RVVs it is 
unknown whether these are truly binary systems, as the radial 
velocity variations may also result from Une profile variability. 
The true value for the observed spectroscopic binary fraction 
in the LEV dataset is thus in the range between 30 ± 6% and 
74 + 6%, depending on how many of the RVVs are true bi- 
nary systems. Note that it is unlikely that the majority of the 
observed RVVs are spurious binaries; the true spectroscopic 
binary fraction is likely close to 74%. 

In this section we recover the orbital size distribution of bi- 
naries in Sco OB2, which can be described with either fa{a) or 
fp(P)- We simultaneously determine the intrinsic binary frac- 
tion Fm- We determine the lower limits flmin and fmin in § 6.1 
and the upper limits Omax and Pmax in § 6.2. The best-fitting 
(power-law) semi-major axis distribution is discussed in § 6.3, 
and the best-fitting (log-normal) period distribution in § 6.4. A 
further discussion on the validity of the independent derivation 
of the orbital size distribution and fq{q) is presented in § 6.5. 
Finally, we summarize our conclusions on the orbital size dis- 
tribution and Fm in § 6.6. 
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Table 7. The five known binaries in Sco 0B2 with an orbital period less than two days. Columns 1-3 list the primary star, the primary spectral 
type, and the measured orbital period. Columns 4 and 5 list the primary mass estimate as derived from the V-band magnitude, under the 
assumption that the mass ratio i?,qx Q and q= \, respectively. Columns 6 and 7 list extremes for the semi-major axis of the binary, as derived 
using Kepler's third law, under the assumption of a mass ratio q ^ Q (6th column), and ^ = 1 (7th column). Finally, colunms 8 and 9 list the 
subgroup of which the binary is a member and the reference for the orbital period. It is possible that a small number of closer, yet undiscovered 
binaries exist in Sco 0B2. 



6. 1. The minimum period and semi-major axis 

The shortest period Pmm can in principle be constrained us- 
ing observations of spectroscopic binaries. Table 7 lists the five 
known binaries in Sco OB2 with an orbital period less than two 
days. These data indicate that the minimum period Pmin is less 
than of order one day. Binary systems with an orbital period 
shorter than one day may be present in Sco OB2, but this or- 
bital period is close to the physical minimum period, leading 
to Roche Lobe overflow. Only a very small fraction of systems 
is expected to have such a short orbital period, as stars in such 
binaries would have been in physical contact during their con- 
traction phase. 

For each short-period binary in Table 7 we have obtained 
an estimate for the semi-major axis, using Kepler's third law. 
We derive the absolute V-band magnitude from the observed 
V-band magnitude using the improved Hipparcos parallax and 
interstellar extinction for each star (de Bruijne 1999). We esti- 
mate a lower and upper mass limit for the primary star, assum- 
ing that the mass ratio is either ^ » or ^ = 1 . As the mass ratio 
of each of these binaries is unknown, we show the correspond- 
ing inferred semi-major axis for the extremes q ^ and ^ = 1. 
As the five short-period binary systems are aU spectroscopic 
binaries, the value of the mass ratio is likely q ^ 0.5. From the 
semi-major axis estimates in Table 7 we thus conclude that Urma 
must be of order 10 R© or smaller. 

As a result of selection effects, the observations cannot be 
used to constrain Prmn and amin any further. However, physical 
arguments can be used to obtain estimates for the lower lim- 
its for f min and amin. The value of amin cannot be significantly 
smaller than the semi-major axis at which Roche lobe overflow 
occurs for one of the components of a binary system. One of 
the components of such a tight binary may fill its Roche Lobe 
if the binary separation a is less than about 2-3 times the radius 
of that star (see, e.g., Eggleton 1983; Hilditch 2001). Assuming 
zero eccentricity, the estimates for amin are of order 12, 4.0, 2.5 
and 1.5 R© for binaries with a primary of spectral type BO V, 
AO V, FO V and GO V, respectively. Adopting a typical mass ra- 
tio of q ^ 0.3, these correspond to minimum orbital periods 
of Pmin = 25, 11, 8, and 6 hours, respectively. In reality the 
latter values are larger due to the presence of eccentric binary 
systems in Sco OB2. A combination of these estimates indi- 
cates that flmin ^ 2 R© and Pmin ^ 12 hours. Binary systems 



with shorter orbital periods are known (see § 2). Although sev- 
eral members of Sco OB2 may have evolved into such a tight 
systems (with P ^ \2 hours), our conclusions on the primor- 
dial binary population wUl not be affected by this assumption. 
If these binaries are present, our inferred binary fraction (§ 8) 
may be mildly underestimated. 

Using observations and physical limitations, we have con- 
strained the values 2 R© 5 amin ^ 10 R© and 0.5 day ^ Pmin ^ 
1 day. Note that these are rather rough constraints. In Section 8 
we will return to this issue when deriving the binary fraction of 
Sco OB2. 

6.2. The maximum semi-major axis and period 

It is difficult to characterize the properties of the widest or- 
bits from observations due to confusion with background stars. 
However, observations of the widest binaries can be used to 
derive a lower limit for the maximum semi-major axis amax. 
Table 8 lists the known binaries in Sco OB2 with angular sep- 
aration larger than 25". Under the assumption that these wide 
binary systems are all physically bound, we additionally list 
a first-order estimate for the semi-major axis a^^i = Dtanp, 
where D is the distance to the binary, and p the angular sepa- 
ration. Table 8 shows that the maximum semi-major axis is of 
order 2 x 10^ R© or larger. Using Kepler's third law we find that 
the orbital period of the binary systems listed in Table 8 is in 
the range between ~ 0.150 Myr and ~ 0.25 Myr (assuming a 
system mass of 5 M© for each binary listed in Table 8). 

A theoretical limit to the maximum semi-major axis Omax is 
obtained from the argument of tidal disruption of binary sys- 
tems. In the Galactic field the maximum observed semi-major 
axis is of the order of ~ 0. 1 pc (5x10^ Rq Bahcafl et al. 1985; 
Close et al. 1990; Chaname & Gould 2004). The analysis of 
Close et al. (1990) shows that ~ 3% of the Galactic disk bina- 
ries has a separation larger than 0.01 pc (0.5 x 10^ R©). Binaries 
in the Galactic halo could be as wide as 1 pc due to the lower 
ambient steUar density (Chaname & Gould 2004). The results 
of the wide binary searches and dynamical simulations above 
indicate that binary systems with a semi-major axis larger than 
a certain value atidai are unstable in the Galactic tidal field and 
are ionized quickly, and that Otidai ~ 0.2 pc (= 9 X 10^ R©) . 

OB associations are expanding groups and are likely un- 
bound (Blaauw 1964; Brown et al. 1999), so that they will 
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Fig. 7. How well does a model with a log-normal period distribution 
fp;ti,a-(P) with centroid // and width cr reproduce the observed angular 
separation distribution fp(fi)l The darkest colors in this figure indi- 
cate the best-fitting models. From dark to light, the contours indicate 
the |cr, Icr, Icr, and 3cr confidence limits for rejection of the model, 
respectively. The cross at (// = 4.8, cr = 2.3) indicates the distribu- 
tion derived by Duquennoy & Mayor (1991) for solar- type stars in the 
solar neighbourhood. The comparison with the observed visual and 
spectroscopic binary fraction is shown in Figure 8. 

dissolve in the field star population within a few tens of Myr. 
If Sco OB 2 is indeed an expanding association, the associa- 
tion must have been denser in the past. The upper limit for the 
semi-major axis may therefore be smaller than atidai- For binary 
systems with a total mass of 5 M©, this tidal limit corresponds 
to ~ 4 Myr. Note that if such a binary system would exist in 
Upper Scorpius, it would have completed only one revolution 
since its birth. For lower-mass binaries, the period correspond- 
ing to a = 0.2 pc would be significantly larger. 

The goal of our study is to find the primordial binary pop- 
ulation, so that theories on star formation can be constrained 
with observations. For our purpose the exact value of amax is 
not a parameter of crucial importance. Binary systems with a 
semi-major axis of order a - 0.2 pc have an orbital period 
larger than ~ 4 Myr. In the context of star formation, the com- 
ponents of these wide binaries may possibly be considered as 
single, as both stars may have formed practically independent 
from each other. 

Combining the information above, we have consttained 2 X 
lO*" Ro ^ flmax ^ 8.9 X 10^ Ro, or alternatively 0.15 Myr <, 
^max ^ 4 Myr. In Section 8 we wiU return to this issue when 
deriving the binary fraction of Sco OB2. 

6.3. A power-law semi-major axis distribution? 

Now that the lower and upper linuts for the semi-major axis 

and period distributions are constrained, it is possible to eval- 
uate these distributions. We first study the possibility that the 
binary population is characterized by a semi-major axis dis- 



tribution of the form /a (a) oc 07", and determine for which 
value of y„ the simulated observations correspond best to the 
observations. We compare the observed angular separation dis- 
tribution fpip) with that of the simulated observations using 
the KS test, and we compare the observed visual and spectro- 
scopic binary fractions with the predictions using the Pearson 
test. In both cases, we test the hypothesis that the observa- 
tions and simulations are realizations of the same underlying 
association model. In our models we adopt semi-major axis 
limits 5 Rq < a < 5x10* Rq, a thermal eccentricity distribution 
fe(e) = 2e and a mass ratio distribution of the form fq(q) oc q'>''i 
with -/q = -0.4. Later in this paper we will discuss the validity 
of these assumptions. 

Figure 6 shows the results of this comparison for models 
with a varying value of y^- The left-hand panel shows the prob- 
ability associated with the KS test when comparing the ob- 
served fpip) and that of the simulated observations. The mod- 
els with -1.15 ^ ja ^ -0.90 produce an observed angular 
separation distribution very similar to that of K05/SHT, and 
consequently have a large goodness-of-fit value. 

The middle panel shows the probability associated with the 
Pearson test when comparing the observed visual binary 
fraction Fm = 31+ 4% and that of the simulated observations 
for models as a function of ya with varying intrinsic binary 
fraction Fm- Models with -0.9 ^ ja ^ -0.6 produce a visual 
binary fraction similar to that of the observations, but only if 
the intrinsic binary fraction is close to 100%. 

The right-hand panel shows the predicted spectroscopic bi- 
nary fraction (SBl, SB2 and RVV) as a function of ja and Fm. 
Note that, unlike the middle panel of Figure 6, we do not show 
the goodness-of-fit between for the observed and predicted 
spectroscopic binary fraction, as we only have lower and up- 
per limits of the former (see below). The two horizontal dotted 
fines in the right-hand panel indicate the lower limit of 30 ± 6% 
(SBl and SB2) and the "upper Umit" of 74 ±6% (SBl, SB2 and 
RVV) for the observed spectroscopic binary fraction Fm- The 
spectroscopic binary fraction increases with decreasing y^, as 
the average orbital separation decreases. Models with ja^ -l 
result in a spectroscopic binary fraction that is compatible with 
the observations. For models with » -1, an intrinsic binary 
fraction very close to 100% is required in order to reproduce 
the lower fimit for the observed spectroscopic binary fraction. 

Considering the results in all three panels of Figure 6, 
we find that the semi-major axis distribution for intermediate 
mass stars in Sco OB2 is best described by /«(«) oc d^" with 
Ja ~ -1.0. The observations are therefore reasonably con- 
sistent with models with Opik's law {ya = -1). Models with 
slightly smaller values of ja correspond somewhat better to 
fpip) and the spectroscopic binary fraction, while those with 
sUghtly larger values of ya correspond better to the observed 
visual binary fraction. The observed visual and spectroscopic 
binary fractions are only reproduced by models with an intrin- 
sic binary fraction close to 100%. 
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Fig. 6. How well can the intermediate-mass binary population in Sco 0B2 be described by a semi-major axis distribution of the form fa(a) oc 
a^" ? Left: the goodness-of-fit for the comparison between the observed angular separation distribution (from the K05 and SHT datasets) and 
that predicted by the models with different values of -/a. A large value for the goodness-of-fit means that the model predictions are consistent 
with the observations. The Icr, 2cr and 3(T confidence limits for rejection of the model are indicated with the horizontal dotted lines. Middle: 
the goodness-of-fit for the comparison between the observed visual binary fraction (from the K05 and SHT datasets) and that predicted by 
the models with an intrinsic binary fraction of 70% (bottom curve), 80%, 90% and 100% (top curve). Right: the predicted (SBl, SB2 and 
RVV) specttoscopic binary fraction (for the LEV dataset) as a function of ja- The five curves indicate an intrinsic binary fraction of 60% 
(bottom curve), 70%, 80%, 90% and 100% (top curve). The bottom horizontal line indicates the lower limit for the observed spectroscopic 
binary fraction, where only the SB Is and SB2s are included, assuming that the radial velocity variations in all RVVs are caused by line-profile 
variability rather than binarity. The top horizontal line indicates the upper limit for the spectroscopic binary fraction, assuming that all RVVs 
are indeed binaries. A combination of the results in these three panels indicates that the semi-major axis disttibution for intermediate mass stars 
in Sco 0B2 can be described by fa{a) oc d'' with exponent ja ~ -1.0 (vertical dotted line in each panel), commonly known as Opik's law. 
An intrinsic binary fraction close to 100% is required in order to produce the large number of visual and spectroscopic binaries detected in 
Sco 0B2. 
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Fig. 8. How well does a log-normal period distribution with centroid /x and width cr reproduce the observed visual and spectroscopic binary 
fractions? Top: the goodness-of-fit for the comparison between the observed and predicted visual binary fraction, as a function of n and cr. Each 

panel represents a set of models with a different intrinsic binary fraction ranging from Fm = 100% (left) to 50% (right). The darkest colors 
indicate the best-fitting models. From dark to light, the three contours indicate the Icr, 2cr, and 3cr confidence limits for rejection of the model, 
respectively. The cross at (/j = 4.8, cr = 2.3) in each panel indicates the distribution derived by Duquennoy & Mayor (1991) for solar-type stars 
in the solar neighbourhood. Bottom: The predicted (SBl, SB2 and RVV) spectroscopic binary fraction as a function of /z, cr and Fm (note that 
these panels do not indicate the goodness-of-fit). The spectroscopic binary fraction is indicated with the gray-shade in each panel, where the 
solid contours indicate the combinations of ji and cr which predict a spectroscopic binary fraction (from black to white) of 80%, 60%, 40% and 
20%, respectively. The other curves indicate the extreme constraints imposed by the observations, of 30 ± 6% (dashed contour) and 74 ± 6% 
(dotted contour), respectively, where the higher value is more likely to represent reality. This figure, combined with the fpip) comparison of 
Figure 7, shows that i/the binary population can be described with a log-normal period disttibution, we require (1) a binary fraction near 100%, 
and (2) a large value for cr (which mimics Opik's law). 
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6.4. A log-normal period distribution? 

We also study the possibility that the binary population can be 
described with a log-normal period distribution (Equation 21) 
for a certain combination of and cr. We simulate associations 
with period distribution fp-^,a-(P) for various values of fi and cr, 
and compare simulated observations with the real observations. 
Again, we consider the observed angular separation distribu- 
tion fpip), the observed visual binary fraction and the observed 
spectroscopic binary fraction. The best-fitting log-normal dis- 
tribution is obtained by combining the results in Figures 7 
and 8, which we describe below. 

In Figure 7 we compare the observed angular separation 
distribution fp(p) for models with varying values of fi and cr. 
The darkest regions in the figure indicate the best-fitting com- 
binations of fi and cr. In general, the models with a large value 
for cr fit the observations well. These models have a very broad 
period distribution, resulting in a /«(«) that is very similar to 
Opik's law, which was shown to be consistent with the ob- 
servations in § 6.3. Models with a period distribution similar 
to that found by Duquennoy & Mayor (1991), i.e., » 4.8 
and cr =s 2.3 are consistent with fpip) as well. Models with 
5 $ /i $ 6 are consistent with the observed distribution fp(p), 
as for these models the peak of logf falls in the visual bi- 
nary regime. This results in an angular separation distribution 
similar to Opik's law. Figure 7 shows that most models with a 
large value of cr are consistent with the observed angular sep- 
aration distribution. These models resemble those with Opik's 
law, which was shown to be consistent with the observations in 
the previous section. Of the models with a large cr, only those 
with extremely wide binary populations produce an in- 

compatible fpip), as only the tail of the distribution then falls 
in the visual regime. 

Figure 8 shows for the same set of models the results for 
the observed visual and spectroscopic binary fractions. The top 
panels show the goodness-of-fit of the visual binary fraction, 
for models with a different value oi fi, a, and for different in- 
trinsic binary fractions Fm (indicated in the bottom-left corner 
of each panel). The darkest regions in each panel indicate the 
best-fitting combinations of fi and cr. The bottom panels show 
the observed spectroscopic binary fraction as a function of 
cr and Fm, where the darkest regions correspond to the highest 
spectroscopic binary fraction. The lower limit (30 + 6%) and 
"upper limit" (74 + 6%) for the spectroscopic binary fraction in 
the LEV sample are indicated with the dotted and dashed con- 
tours, respectively. Any well-fitting model has a combination 
of jj. and cr within these contours, preferably near the left-hand 
contour, as a significant number of the RVVs detected by LEV 
is likely to be a true binary. In each panel of Figure 8 we in- 
dicate the values (/i, cr) - (4.8,2.3) found by Duquennoy & 
Mayor (1991) for solar-type stars in the solar neighbourhood. 

By comparing the top and bottom panels in Figure 8 we find 
that the intrinsic binary fraction is larger than « 70%, whatever 
the values of /i and cr are. Models with Fm ~ 100% are most 
consistent with the observations. The Duquennoy & Mayor 
(1991) period distribution with (yU,cr) = (4.8,2.3) is consis- 
tent with the observed angular separation distribution and vi- 
sual binary fraction, but underpredicts the number of spectro- 



scopic binaries, even if the intrinsic binary fraction is 100%, 
and all RVVs detected by LEV are spurious. Similar models 
with a larger value for cr are consistent with the observations. 
However, these period distributions are very broad, making the 
binary population indistinguishable from that resulting from 
Opik's law. 

6.5. The ambiguity in deriving the orbital size 
distribution and fq(q) 

In Sections 5 and 6 we have constrained the orbital size dis- 
tribution and the mass ratio distribution fgiq) independently. A 
risk associated with this approach is that one may find a local, 
rather than a global set of best-fitting solutions. In this section 
we demonstrate that our choice has been appropriate, and that 
we have indeed found the best-fitting model. 

Even if the orbital size distribution and fgiq) are intrinsi- 
cally independent, there could be a correlation between a (or P) 
and q among the detected binaries. The latter is important for 
the analysis of the results of a spectroscopic binary survey (e.g., 
Kobulnicky et al. 2006), but is of less importance for visual bi- 
naries. In a visual binary survey the orbital size is related to the 
angular separation p, and q is related to the brightness contrast 
between primary and companion. As the maximum detectable 
magnitude difference depends on the separation between the 
binary components, a correlation between the orbital size and 
q may be introduced among the closest detected companions. 

In order to study whether we have indeed found the global 
set of best-fitting solutions, we extend the parameter space, and 
make comparisons for models with varying orbital size dis- 
tribution and varying fqiq). We hmit ourselves to mass ratio 
distributions of the form fqiq) cc q^i . We compare the two- 
dimensional distribution fp,qip, q) and the visual binary fraction 
Fm of the simulated observations of the combined K05/SHT 
dataset with that of the observations. We compare the observed 
distribution fp,qip, q) using the two-dimensional KS test (e.g., 
Press et al. 1992), and we compare the observed visual binary 
fraction using the Pearson test. 

Figures 9 and 10 show plots corresponding to those in 
Figures 5 and 6, but now with the additional free parameter 
jq. The results for a power-law semi-major axis distribution 
in Figure 9 are similar to those derived in Sections 5 and 6. 
The best-fitting values for jg and jq in these diagrams are 
equal (within the error bars) to those of the independent deriva- 
tions of faia) and fqiq). Note that the adopted combination 
ija ~ -1-0,74 ~ -0.4) does not fall in the darkest regions 
of Figure 9. The comparison for fp,qip, q) (left panel) suggests 
somewhat smaller values of ja and y^, while the comparison 
for the visual binary fraction suggests larger values. An inspec- 
tion of both comparisons suggests that our adopted combina- 
tion (7a * - 1 .0, Jq ~ -0.4) is among the best-fitting combina- 
tions for both fp,qip,q) and Fm. The figure also clearly shows 
again that an intrinsic binary fraction close to 100% is required. 
We do not attempt to further constrain the combination (7,,, jq), 
as this would involve a combination of the independent com- 
parisons for fp,qip,q) and Fm, as this would require accurate 
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Fig. 9. How well do models with a semi-major axis distribution fa{a) oc a'" and a mass ratio distribution fqiq) tx q'''' correspond to the 
observations? This figure is the multi-dimensional equivalent of Figures 5 and 6, and indicates that with the independent derivation of 
and ja we have found a globally best-fitting solution. Left: the consistency between the observed two-dimensional distribution fp^ip, q) and 
that of simulated observations, for models with different values of -y„ (horizontal cais) and jq {vertical axis), respectively. The darkest colors 
indicate the best-fitting models. From dark to light, the three contours indicate the Icr, 2cr, and 3cr confidence limits for rejection of the model, 
respectively. The cross in each panel indicates the values of and y, that we adopt to describe the properties of Sco OB2. Right: the consistency 
between the observed visual binary fraction and that of the simulated observations, for models with different values of y^ and y^. Each panel 
corresponds to a set of models with a different binary fraction, which is indicated in the bottom-left comer. The darkest colors indicate the 
best-fitting models. From dark to light, the three contours indicate the Icr, 2cr, and 3cr confidence limits for rejection of the model, respectively. 
Models with a binary fraction smaller than » 100% are inconsistent with the observations, as they underpredict the observed number of visual 
binaries. The cross in each panel indicates the best-fitting values of ya « -1.0 and y, » -0.4 for Sco 0B2. 



knowledge of the relative weights that should be assigned to 
each of these comparisons. 

The results for the log-normal period distribution /p;^,o-(P) 
shown in Figure 10 also indicate that our independent deriva- 
tion of fp-^,a-{P) and fqiq) is justified. The top panels show that 
the models with a mass ratio distribution exponent = -0.7, 
-0.6, -0.3 and -0. 1 produce a distribution fp,q{p, q) that is less 
consistent than that for the models with the best-fitting solution 
jq = -0.4, for any combination of /i and cr. The bottom panels 
of Figure 10 indicate for the observed visual binary fraction the 
best- fitting combinations of ju, cr, jq and Fm (darkest regions). 
These comparisons indicate that, if the log-normal period dis- 
tribution with (ju, cr) - (4.8, 2.3) holds, the intrinsic binary frac- 
tion should be Fm ^ 70%. Again, we do not attempt to find the 
optimum combination (//, cr, jq, Fm) by combining the results 
of all panels, as this would require a detailed prescription of 
the relative weights given to the comparisons of fp,q{p, q) and 
the visual binary fraction, which is not trivial. Note that in our 
models we adopt a binary fraction that is independent of semi- 
major axis and mass ratio. 

6.6. Conclusions on fa(a), fp(P) and Fm 

In the sections above we have constrained the properties of the 
orbital size distribution among intermediate-mass binaries in 
Sco OB2, using the visual binary surveys of K05 and SHT and 
the spectroscopic binary survey of LEV. The tightest binaries 
have IRq ^ amm ^ 10 R©, corresponding to 0.5 day ^ Pmm ^ 



1 day. The widest orbits have 2 x lO*" R© ^ amax ^ 8.9 x 
10^ Ro, or 0.15 Myr < Pmax ^ 4 Myr. We have compared the 
observations with models that have either a power-law semi- 
major axis distribution or a log-normal period distribution, in 
order to find the best-fitting orbital size distribution. 

We considered a power-law semi-major axis distribution 
faifl) with varying exponent ja- Models with a semi-major axis 
distribution exponent ja ~ -1.0 are consistent with the ob- 
served angular separation distribution and the observed visual 
binary fraction, but only if the binary fraction is close to 100%. 
We have additionally considered models with a log-normal pe- 
riod distribution fp-^i^a-iP) with varying centroid // and width tr. 
The distribution derived by Duquennoy & Mayor (1991), with 
H = 4.8 and cr - 2.3 is shown to be consistent with the ob- 
served angular separation distribution and the observed visual 
binary fraction for models with an intrinsic binary fraction near 
100%. The observed spectroscopic binary fraction, however, is 
rather large compared to those predicted by the above models, 
even if all stars are in binary systems. 

Opik's law could be considered as resulting from a very 
broad log-normal period distribution (cr oo). Although 
Opik's law and the log-normal period distribution with (fj., cr) = 
(4.8, 2.3) are consistent with the observations (i.e., they cannot 
be excluded), the observations indicate that the best descrip- 
tion of the binary population in Sco OB2 is intermediate, with 

2.3 < cr < oo. 

We have constrained the mass ratio distribution and the or- 
bital size distribution separately. Our investigation of this ap- 
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Fig. 10. How well do models with a log-normal period distribution fp-^,a-(P) and a mass ratio distribution fg(q) oc q'l correspond to the 
observations? Top: a comparison for the observed two-dimensional distribution fp,q(p, q) for models with different values of jx, cr and y^. The 
cross in each panel indicates the values observed by Duquennoy & Mayor (1991): ji = 4.8 and cr = 2.3. The four panels correspond to models 
with a mass ratio distribution exponent = -0.7 {left), -0.6, -0.3 and -0.1 (right), respectively (cf. Figure 8). The goodness-of-fit is best for 
the dark-shaded regions. Contours are plotted at Icr, 2cr and 3cr confidence regions for rejection of the model. Bottom: the comparison for the 
observed visual binary fraction for models with different values of yu, cr and y^. The binary fraction, indicated in the bottom-left corner of each 
panel, decreases from 100% (left) to 50% (right). Contours are plotted at Icr, 2cr and 3cr confidence regions for rejection of the model. 



proach in section 6.5 indicates that we have indeed found the 
global best-fitting model. 

7. Recovering the eccentricity distribution 

The eccentricity distribution for binaries in Sco OB2 is not 
known. The field star population, on the other hand, has been 
studied in detail by Duquennoy & Mayor (1991), who sur- 
veyed solar-type stars in the solar neighbourhood for bina- 
rity and find the following composite eccentricity distribution. 
The close binaries (P < 10 days) have a negligible eccen- 
tricity, presumably due to tidal circularization, although they 
may have formed in nearly circular orbits. Binaries with a pe- 



riod 10 < P < 1000 days show a bell-shaped eccentricity 
distribution that peaks near e - 0.3. Binaries with a larger 
orbital period tend to have a thermal eccentricity distribution 
fe(e) = 2e, though the presence of large-eccentricity binaries 
are inferred by incompleteness corrections. The latter distribu- 
tion is expected from energy equipartition (Heggie 1975) as a 
result of multiple soft encounters. A similar result was found 
by Mathieu (1994), who studied the properties of pre-main- 
sequence binaries. 

Both studies indicate that the eccentricity distribution may 
well be a function of orbital period, which is perhaps also the 
case for Sco OB2. The very tight binaries tend to have more 
or less circular orbits. Tidal circularization may have affected 
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these binaries in Sco OB2. The tidal circularization period Pciro 
defined as the orbital period at which a binary orbit with a typ- 
ical initial orbital eccentricity circularizes, is Pcirc = 7-10 
days for young stellar populations (Meibom & Mathieu 2005). 
The intermediate-period binaries have a somewhat flat eccen- 
tricity distribution, possibly reflecting the formation process. 
The widest binaries tend to have a more thermal eccentricity 
distribution, which may (partially) be a result of binary-binary 
interactions. It must be noted, however, that the eccentricity 
distribution for the latter two period ranges is very difficult to 
measure because of the large orbital periods. 

It is practically impossible to derive the properties of the 
eccentricity distribution from observations of visual binaries 
(e.g., Kouwenhoven 2006). On the other hand, this means that 
the results we derived from the observations of visual binaries 
(i.e. fq{q), fa(a), fp(P) and Fm) are unaffectedhy our choice for 
/e(e) in the models. In principle it is possible to constrain fe(e) 
using observations of spectroscopic and astrometric binaries. 
However, only one Sco OB2 binary with an astrometricaUy de- 
termined orbital solution is known. 

Only the measurements of the spectroscopic binaries can 
thus be used to constrain the eccentricity distribution for bina- 
ries in Sco OB2. The properties of the eccentricity distribution 
derived below are only valid for the short-period (P ^ 14 days; 
a ^ 40 R©) and intermediate-mass binaries. Caution should be 
taken with generalizing the results for the full population of 
Sco OB2. As we do not model the detection of SB Is, SB2s and 
RVVs in detail, we cannot compare the eccentricity distribution 
resulting from the simulated observations (SB Is, SB2s, and 
RVVs) with the observed eccentricity distribution (SB Is and 
SB2s only). However, it is possible to set several constraints 
on the intrinsic eccentricity distribution using the observations 
of LEV. Due to the Umited number of eccentricity measure- 
ments, the complicated selection effects, and the possible cor- 
relation between e and P, we do not attempt to fully recover 
feie). Instead, we consider three commonly adopted eccentric- 
ity distributions, and compare the resulting predictions with the 
observations. 

In our analysis we will consider three frequently used 
eccentricity distributions: the flat eccentricity distribution 
/flat(e) = 1, the thermal eccentricity distribution f2e(e) - 2e, 
and the single- valued eccentricity distribution fe„(e) - 6(e-eo) 
with < Co < 1- Simply by analyzing Figure 2, we can 
rule out an eccentricity distribution of the form fe^ , where all 
binaries have the same eccentricity eo. For associations with 
feie) = /co(e) the distribution fe(e) is unbiased, as all binaries 
have eg- An error is associated with each eccentricity measure- 
ment, so that fe(e) is broader than fde). In the spectroscopic 
binary sample, the error in the eccentricity is of order 0.05, rul- 
ing out the best-fitting distribution with eo = 0.27 with ~ Bo- 
confidence. 

The observed eccentricity distribution is in better agree- 
ment with the thermal distribution fze and the flat distribution 
/flat- For these distributions the relative (intrinsic) fraction of 
binary systems with e < 0.5 is: 



„ # binaries with e < 0.5 

Ele = TTT-. ■■ = 25% , (23) 

# bmaries 



En-it - 



# binaries with e < 0.5 



50% 



(24) 



# binaries 

In the LEV sample, 14 of the 16 targets for which the orbital 
elements are determined have e < 0.5, and only 2 have e > 0.5. 
The apparent overabundance of low-eccentricity (e < 0.5) sys- 
tems (see also Figure 2) can partially be explained by selection 
efl'ects. Highly eccentric systems spend a large fraction of their 
orbit near apastron, and are therefore more difficult to detect. 

In § 6 we have shown that the observed distribution fp(p) 
and binary fraction Fm of visual binaries are consistent with 
both Opik's law and a log-normal period distribution with 
Fm ~ 100%. For each of these distributions we predict the 
number of spectroscopic binaries (SB Is, SB2s and RVVs) in 
the LEV dataset. For models with Opik's law, our simulations 
indicate that LEV would be able to detect ~ 38% of the bi- 
nary systems spectroscopically (i.e., as SBl, SB2 or RVV) for 
a flat eccentricity distribution, and ~ 36% for a thermal ec- 
centricity distribution. If the log-normal period distribution of 
Duquennoy & Mayor (1991) holds, we find ~ 27% for a flat 
eccentricity distribution, and ~ 24% for a thermal eccentricity 
distribution, respectively. The observed spectroscopic binary 
fraction is thus only mildly dependent on fde), and decreases 
with an increasing average orbital eccentricity. The spectro- 
scopic binary fraction is also proportional to Fm- The expected 
fraction of spectroscopic binary systems with e < 0.5 among 
the observed targets is thus given by: 



EiejDM = 24% X Fm X < 6.0% ; 

£2.,6pik =36%xFmx£2. <9.0%; 

Fflat,DM = 27% X Fm X Ea^t < 13.5% ; 

Fflat,op.k = 38% X Fm X Ea,t < 19.0% , 



(25) 
(26) 
(27) 
(28) 



where the upper limits on the right-hand side are obtained for 
models with an intrinsic binary fraction Fm = 100%. 

For each of the 23 RVVs in the LEV sample we do not 
know whether it is truly a binary, and if so, we do not know its 
eccentricity. We therefore constrain f^ie) by considering two 
extreme cases for Flev- (1) none of the RVVs have e < 0.5 
(i.e., they are either spectroscopically single, or they are bina- 
ries with e > 0.5), and (2) all RVVs are binary systems with 
e < 0.5. In these extreme cases, Flev is constrained by 



26 ± 6% < Flev < 70 ± 6% . 



(29) 



The observed spectroscopic binary fraction Fm is inconsistent 
with the predictions in Equations 25-28. The models with the 
log-normal period distribution (with = 4.8 and a = 2.3) or 
the thermal eccentricity distribution do not reproduce the obser- 
vations; they underpredict the number of binaries in the LEV 
sample with e < 0.5, even if all RVVs detected by LEV are spu- 
rious (which is very unlikely), and the model binary fraction is 
100%. The model with Opik's law, a flat eccentricity distribu- 
tion, and a binary fraction of 100% (Equation 28) is consistent 
with the observations, but only if most of the RVV candidate 
binaries are (spectroscopically) single stars exhibiting line pro- 
file variabiUty. 

The analysis above indicates that the short-period, 
intermediate-mass binaries do not follow the thermal eccentric- 
ity distribution, but have lower eccentricities on average (see. 
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however, § 10.1). Note that the discussion on the eccentricity 
distribution is constrained only by a small number of short- 
period binaries with a massive primary. In the LEV sample, 
12 out of the 16 systems with a measured eccentricity have a 
period P < Pdrc, suggesting that the derived properties of the 
eccentricity distribution may not be representative for the bi- 
nary population of Sco OB2. Accurate radial velocity and as- 
trometric surveys among intermediate- and low-mass members 
of Sco OB2 are necessary to further characterize the eccentric- 
ity distribution. 

8. The binary fraction in Sco 0B2 

In the previous sections we have constrained the pairing func- 
tion, the mass ratio distribution, the orbital size distribution 
and the eccentricity distribution. The analysis of each of these 
properties indicates that the binary fraction among intermediate 
mass stars in Sco OB2 has to be close to 100%. 

Table 9 lists the fraction of detected binaries among the tar- 
geted sample, for the different surveys. Results are listed for the 
observations and for the simulated observations of models with 
Opik's law and for those with the Duquennoy & Mayor (1991) 
period distribution. For each of these models we have adopted 
an intrinsic binary fraction Fu - 100%. The values listed in 
Table 9 are upper hmits, and are proportional to Fm- 

The observed visual binary fraction is consistent with 
Opik's law, but also with the log-normal period distribution 
(ju = 4.8, cr = 2.3). A binary population with intermediate 
properties (2.3 < cr < oo), however, is most consistent with 
the observations. The observed spectroscopic binary fraction 
cannot be compared directly with the simulated observations, 
as the true nature of the detected radial velocity variables is 
unknown. The observed spectroscopic binary fraction can be 
Umited by 30 ± 6% (in the unUkely case that all RVVs are spu- 
rious) and 74 ±6% (if all RVVs are indeed binaries). Even if the 
intrinsic binary fraction is 100%, models with Opik's law only 
marginally reproduce the observed spectroscopic binary frac- 
tion, while those with the Duquennoy & Mayor (1991) period 
distribution underpredict the observed value. In the case that 
50% or 100% of the RVVs are indeed binaries, the log-normal 
period distribution can be excluded with 4cr and 5cr confidence, 
respectively, if om other assumptions hold (see § 10.1). 

Our model for the Hipparcos selection effects for the (X), 
(O), and (G) binaries is not strict enough. This is because in 
our models we consider all binaries with the properties listed in 
Table 6 as astrometric binaries. In reahty, however, only a sub- 
set of these would have been marked as an astrometric binary 
by Hipparcos (see § 4.6.2). The predicted astrometric binary 
fraction can therefore only be used as an upper Umit when a 
comparison with the observations is made. The predicted frac- 
tion of visual binaries in category (C) on the other hand, is well- 
modeled. Models with a binary fraction Fm ~ 100% correctly 
predict the number of (C) binaries. 

An inspection of the observed binary fractions for the dif- 
ferent datasets (Table 9) and their associated statistical errors 
above indicates that the binary fraction among intermediate 
mass stars in Sco OB2 must be larger than ^ 85% at the 2cr 
confidence level, and larger than « 70% at the 3<t confidence 



level. We find that models with a binary fraction of 100% are 
most consistent with the observations; models with a signif- 
icantly smaller binary fraction are inconsistent with the ob- 
served binary fraction. Mason et al. (1998) found that almost 
all high-mass stars are multiple. Our complementary results in- 
dicate that also practically all intermediate mass stars form in a 
binary or multiple system. 



8.1. Dependence of the derived binary fraction on the 
limits of a and P 

For the models with a power-law semi-major axis distribu- 
tion faia) we have adopted a semi-major axis distribution 
5Ro < a < 5x1 0^ Ro, and for models with a log-normal 
period distribution fp-^i^^a-iP) we have adopted the period range 
0.7 day < P < 3 x 10^ day. These lower and upper limits 
are derived from the observational data and have their associ- 
ated uncertainties (see § 6.1 and 6.2). In this section we briefly 
study the efl'ect of these uncertainties on our results. 

The binary population properties derived from the angu- 
lar separation distribution of visual binaries are not affected by 
our choice of the lower and upper limits. The results in the left- 
hand panels in Figures 5, 6, 8, 9, the top panels in Figure 10, 
and the derived distributions /«(«), fp(P) and fqiq) are insensi- 
tive to our choice for the tightest or widest orbits. The inferred 
binary fraction, however, does depend on these limits. Table 10 
shows the expected visual and spectroscopic binary fractions, 
for models with different Omin and amax (adopting Opik's law), 
and different Pmm and Fmax (adopting the period distribution of 
Duquennoy & Mayor 1991). If, for a given model, we increase 
the value of amax or Pmax. this leads to smaller values for the 
simulated visual and spectroscopic binary fractions. If we in- 
crease flniin or f min. this Icads to Smaller spectroscopic binary 
and a larger visual binary fraction. In other words, an overes- 
timation of amax or f max leads to an overestimation of the in- 
ferred intrinsic binary fraction Fu- An overestimation of Omin 
or -T min leads to an overestimation of Fm as derived from the 
observed spectroscopic binary fraction, but to an underestima- 
tion of Fm as derived from the observed visual binary fraction. 

Fortunately, the properties of the tightest and widest bina- 
ries are reasonably well constrained by observations. Table 10 
shows that the uncertainty in the observed visual binary frac- 
tion is approximately 2%, as a result of the errors in amin and 
flmax- The corresponding error in the spectroscopic binary frac- 
tion ranges from about 2% for the log-normal period distribu- 
tion to approximately 5% for Opik's law. The uncertainty in the 
inferred intrinsic binary fraction is slightly smaller as we com- 
bine the visual and spectroscopic results. An overestimation of 
«min, for example, leads to a larger visual binary fraction but 
a smaller spectroscopic binary fraction, so that the systematic 
uncertainty on Fm partially cancels out. The effects of the un- 
certainty in the lower and upper hmits of fa{a) and fp{P) are 
therefore limited. We estimate that this uncertainty results in a 
systematic error of ;S 2 - 4% in the inferred intrinsic binary 
fraction Fm. 
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Dataset 


/ NU 

(JDS. 


fraiaXya = -1 






= 4.0, cr = 


z.i 








Fuilo-) 


FM(2cr) 


FM(3cr) 


Pm 


FmCIct) 


Fm(2o-) 




K05 


30 ±4 


26 


« 100 


91 - 100 


81 - 100 


35 


97 - 100 


70 - 100 


61 - 100 


K06 


< 82 ± 8 


27 








37 








SHT 


26 + 4 


23 


94 - 100 


77 - 100 


63 - 100 


34 


64-97 


52 - 100 


43 - 100 


LEVSBl/SB2only 


> 30±6 


36 


> 68 


>54 


>42 


24 


« 100 


> 81 


>64 


LEVSB1/SB2/RVV 


< 74± 6 


36 








24 








BRV analysis - RVV 


< 60±5 


39 








27 








BRVSB1/SB2/RVV 


<66±5 


39 








27 








HIP (X)/(0)/(G) 


4± 1 


< 31 


> 12 


> 10 


> 8 


<25 


> 14 


> 11 


> 9 


HIP (X)/(0)/(G)/(S) 


9+ 1 


< 31 


> 25 


> 22 


> 19 


< 25 


> 30 


> 26 


> 22 


HIP(C) 


15 ± 2 


13 


- 100 


96 - 100 


86- 100 


18 


77 - 94 


68 - 100 


61 - 100 



Table 9. The observed binary fraction (in %) and inferred intrinsic binary fraction (in %) for the different datasets discussed in this paper. 
Columns 1 and 2 list the dataset and the observed binary fraction. The predicted observed binary fraction resulting from Opik's law for 
each dataset (adopting Fm = 100%) is listed in column 3, followed by the inferred Icr, 2cr and 3cr confidence ranges of the inferred binary 
fraction. The predicted observed binary fraction for the log-normal period distribution with fi = 4.8 and cr = 2.3 (adopting Fm = 100%) 
and corresponding confidence ranges for the intrinsic binary fraction are listed in columns 7-10. The adopted association parameters are 
listed in Table 1. For models in columns 3-6 the semi-major axis range is 5 Rq < a < 5 x 10* Rg. For models in columns 7-10 the period 
range is 0.7 day < P < 3 x 10* day. For each model we assume a mass ratio distribution of the form fg{q) oc with = -0.4 and a 
thermal eccentricity distribution. The comparison between observations and simulated observations indicates that the binary fraction among 
intermediate-mass binaries in Sco 0B2 is close to 100% (^ 70% with 3cr confidence). 



/opik(«) hmits 


Fm,K05,SHT (%) 


Fm.lev (%) 


2 Ro 2.0 X 10* Ro 


25 


42 


2 Ro 8.9 X 10« Ro 


24 


37 


10 Ro 2.0 X 10* Ro 


27 


34 


10 Ro 8.9 X 10* Ro 


25 


30 


/dm(P) limits 


Fm VB (%) 


Fm SB (%) 


0.5 day 0.15 Myr 


38 


25 


0.5 day 4.00 Myr 


35 


22 


2.0 day 0.15 Myr 


39 


22 


2.0 day 4.00 Myr 


36 


20 



Table 10. The simulated observed binary fraction for models with 
a different assumption for the tightest and widest binaries. The first 
and second column list the properties of the tightest and widest or- 
bits, respectively. The third and fourth colunm list the visual binary 
fraction for the simulated combined K05/SHT observations and the 
spectroscopic binary fraction for the simulated LEV observations, 
respectively. For each model we adopt an intrinsic binary fraction 
Fm = 100% and a thermal eccentricity distribution. The values listed 
in columns 3 and 4 are proportional to Fm. The statistical errors on 
the listed binary fraction are 0.5 - 1%. The values listed in this table 
provide error estimates for our results in Table 9. The uncertainties in 
the Umits of fa(a) and fp{P) result in a systematic error of < 2 - 4% 
in the inferred intrinsic binary of the binary population. 



9. The primordial binary population in Sco 0B2 

Sco OB2 is a young OB association (5-20 Myr) with a low stel- 
lar density (~ 0.1 Mopc"-'), comparable to the stellar density 
in the solar neighbourhood. One therefore expects that stellar 
and dynamical evolution have only mildly altered the proper- 
ties of the binary population. In this section we study which 
binary systems in Sco OB2 may have changed one or more of 
their parameters since the time of formation. We first consider 
whether stellar (and binary) evolution has affected the binary 
population, and subsequently investigate the importance of dy- 
namical evolution. 



For stellar populations with an age less than about 20 Myr 
the fraction of binaries of which the properties have changed 
due to stellar evolution is small. The hydrogen-burning time 
as a function of initial mass and initial rotational velocity was 
studied by Meynet & Maeder (2000). According to their model, 
stars more massive than 32 - 42 M© (06V) in US, and stars 
more massive than 10 - 11 M© (B2V) in UCL and LCC have 
evolved away from the main sequence. The lower and upper 
limits for each subgroup correspond to stars with an initial ro- 
tational velocity of kms"' and 300 kms"', respectively. The 
most massive star in the US subgroup is Antares, an Ml lb su- 
pergiant with an initial mass of 22.5 ± 2.5 M© (Preibisch et al. 
2002), indicating that we may have mildly overestimated the 
upper mass limit for the US subgroup. The derived mass limit 
for the UCL and LCC subgroups corresponds roughly to the 
observed range of spectral types. An extrapolation of the mass 
distribution (Equation 2) to infinity suggests that only 3-8 stars 
of spectral type O (M ^ 20 M©) have formed in Sco OB2. 

Several of the binary systems with a primary initially more 
massive than the tum-ofl" mass will have produced a compact 
component. These compact objects and their previous compan- 
ion stars may have obtained a significant kick during the super- 
nova event, resulting in a space velocity significantly different 
from that of their parent association (runaway stars). The prop- 
erties of compact binary systems including at least one star with 
a higher initial mass may have changed due to binary evolution. 
For ^ 99.9% of the binaries in Sco OB2, however, both com- 
ponents have an initial mass (much) less than those mentioned 
above; stellar evolution will not have afl'ected these binaries. 

Changes of the binary population due to stellar evolution 
can therefore be neglected. The change due to dynamical evo- 
lution is more complicated to quantify. Below we calculate 
which binary systems are afi'ected by dynamical evolution, at 
the current stellar density of the association. After that we 
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briefly discuss the consequences of the fact that the subgroups 
of Sco OB2 may have been denser in the past. 

The subgroups of Sco OB2 can be approximated as roughly 
spherical, with radii of order /? = 20 - 30 pc (see Table 1). 
The structure and kinematics of Sco OB2 have been studied 
by de Bruijne (1999), who finds that the one-dimensional in- 
ternal velocity dispersion for each of the subgroups is (Tv,id ^ 
1.0 - 1.5 kms"', corresponding to a three-dimensional veloc- 
ity dispersion 0-^30 ^ 1.7 - 2.6 kms ^ The crossing time of 
each subgroup is thus Tcross = ^/Cv.sd ~ 10 Myr. Preibisch 
et al. (2002) find that the US subgroup contains approximately 
2525 members more massive than 0.1 M©. The total number 
of members A'^ in the substellar regime depends on the un- 
known values of Mmin and Mp in Equation 2. For reasonable 
value of Mnun and Mp, about 44% of the members have a 
mass smaller than 0.1 M©, so that the total number of sys- 
tems equals ^ 4500. Adopting a mass ratio distribution of 
the form fq{q) oc ^ and a binary fraction of 100%, the 
total mass of US is of order 2280 Mq. We assume that the 
UCL and LCC subgroups contain a similar number of mem- 
bers. The mass density of each subgroup is therefore of order 
0.04 ± 0.02 Mq pc"^, corresponding to 0.09 ± 0.04 systems (or 
0.17 ± 0.09 individual stars) per cubic parsec. 

An estimate for the relaxation time Treiax> i-e., the typical 
timescale at which a system of A'^ single stars has "forgotten" 
the initial conditions, is 



'^relax 



81nA^ 



^ 670 Myr 



(30) 



(Binney & Tremaine 1987). As the cross-section for binary sys- 
tems is significantly larger than that of single stars, they interact 
more frequently, so that in reality the timescale of dynamical 
evolution is shorter than given by Equation 30. Whether a bi- 
nary system experiences a strong encounter, depends not only 
on the properties of the association, but also on the size (more 
specifically, the binding energy) of the binary systems. A bi- 
nary system is often classified as hard or soft, depending on 
whether it is likely to experience a strong encounter within a 
relaxation time. A binary system is called "hard" if its orbital 
energy Eorb is significantly larger than the mean kinetic energy 
{K) of the surrounding stars: fiorb > 3(/r) (Heggie 1975; Hills 
1975); otherwise a binary system is called "soft". The orbital 
energy of a binary system is given by 



^orb — 



qGM\ 
la 



(31) 



where q is the mass ratio, G the gravitational constant. Mi the 
primary mass, and a the semi-major axis of the binary system. 
The mean kinetic energy is given by 



(32) 



where {Mj) is the average mass of a field object, and iTi^d its 
three-dimensional velocity dispersion. Assuming a binary frac- 
tion of 100%, the median mass of a field object (i.e., another bi- 
nary system in the association), is {Mj) « 0.21 Mq. The three- 
dimensional velocity dispersion is 0-^30 ^ 1 .7 - 2.6 km s ' (de 
Bruijne 1999). The hard/soft boundary ^hs for a binary star with 



primary mass Mi and mass ratio q is thus given by ^orb = 3{X'), 
i.e.. 



aha = 



■qM\ 



v,3D 



(33) 



^350^(^) AU = 7.5xl04^(^) R©. 



The latter value is in good agreement with the widest known bi- 
naries in Sco OB2 (see § 6.2), which all have a primary mass of 
order 3 M©. The datasets studied in this paper contain mainly 
A and B type members of Sco OB2, which have Mi S 1.4 M©. 
Among the binaries more massive than 1.4 M©, about 72 ± 5% 
are hard, while the widest 28 + 5% of the binaries are soft, 
if Opik's law holds. If the Duquennoy & Mayor (1991) pe- 
riod distribution holds, 77 ± 3% of these binaries are hard, and 
24 ± 3% are soft. The quoted errors include only the contri- 
bution from the uncertainty in a^^^ and flmax- The population 
of A and B stars must therefore be close to primordial, given 
the young age of the association. On the other hand, among the 
lower-mass binaries (Mi « 0.5 M©) about half of the binaries is 
hard. The other half of the binaries (i.e., the soft low-mass bina- 
ries) is expected to experience a close encounter within a relax- 
ation time. However, as shown above, the (current) relaxation 
time is significantly larger than the age of Sco OB2, suggesting 
that this population is unlikely to have changed significantly. 
The majority of binaries, in particular the intermediate-mass 
binaries studied in this paper, are unhkely to have experienced 
a strong encounter during the lifetime of the association, as- 
suming that the density did not change over time. 

The density of Sco OB2 may have been (much) higher dur- 
ing the star formation process. Simulations of van den Berk 
et al. (2007), however, suggest that even at higher stellar den- 
sities, the binarity population is only mildly afl'ected by dy- 
namical interactions. In the latter paper, simulations of small 
{N = 100) and initially dense (~ 10^ systems per cubic parsec) 
expanding star clusters are presented. Their simulations include 
triple systems for which the period of the outer orbit reaches up 
to 1000 years. Although these simulated clusters are more than 
a milUon times denser than the current density of Sco OB2, 
the properties of the binary population do not change signifi- 
cantly within 20 Myr. As about 60% of the binary systems in 
Sco OB2 has an orbital period smaller than 1000 years, these 
binaries are expected not to have changed significantly since 
the birth of the association, even if the density of Sco OB2 was 
orders of magnitude larger at the time of formation. 

The discussion above indicates that the current binary pop- 
ulation is likely very similar to the current binary population, 
suggesting that all intermediate mass stars have formed in a bi- 
nary or multiple system. However, due to the as yet unknown 
initial conditions of Sco OB2 (such as the initial stellar den- 
sity), N-body simulations of expanding OB associations are 
necessary to further constrain the primordial binary population 
of Sco OB2. By varying the initial conditions, evolving each 
simulated association for 5-20 Myr, and comparing the out- 
come with the observations, the primordial binary population 
can be recovered. The latter technique is referred to as inverse 
dynamical population synthesis (see Kroupa 1995a,b). 
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10. Discussion 

10.1. Triple and higher-order systems 

In our analysis we have constrained the properties of the bi- 
nary population primarily using observations of visual bina- 
ries. Our simulations have indicated that even models with an 
intrinsic binary fraction of 100% produce a rather low spec- 
troscopic binary fraction as compared to the observations of 
LEV. Throughout this paper, however, we have ignored the 
presence of triple and higher-order systems. These systems are 
known to be present in Sco OB2: the observed higher-order 
multiplicity fraction among Hipparcos members of Sco OB2 is 
{T + Q + . ..)/(S + B + T + Q + ...) = 8.4% (see Table 1). 
This value is a lower limit due to the presence of undetected 
companions. The studies of Tokovinin & Smekhov (2002) and 
Correia et al. (2006) have indicated that 20 - 30% of the wide 
visual binaries have a spectroscopic subsystem. 

The presence of these triple and higher-order systems 
among the Hipparcos members of Sco OB2 could explain the 
apparent underabundance of spectroscopic binaries in our mod- 
els at least partially, if one would only include the outer com- 
ponents of a multiple system in the statistics. We do not attempt 
to calculate the contribution of spectroscopic subsystems in this 
paper, as this would require detailed a-priori knowledge about 
the triple population. Further detailed observational studies are 
necessary to characterize the properties of these systems, and 
to derive the primordial binary-and-multiple-systems popula- 
tion of Sco OB2. 

10.2. Comparison with Heacox (1995) 

Heacox (1995) derived the mass ratio distribution for Sco OB2 
binaries with an intermediate-mass primary from the LEV 
dataset. He derived the mass ratio distribution using the ob- 
servations of all 22 binaries with spectroscopic elements, re- 
sulting in a mean mass ratio (q) « 0.27 with a standard deriva- 
tion cTg » 0.04. In our analysis we adopt a mass ratio distri- 
bution of the form fg(q) oc q^'i. If this distribution is adopted, 
the mean mass ratio is (q) = {y^j + l)/(yq + 2), and the stan- 
dard deviation is cr^ = (y^ + l)/iyq + 3) - (q)^. Our best-fitting 
solution has yg = -0.4 ± 0.1, so that (q) = 0.37 ± 0.04 and 
CTq = 0.090 + 0.002 (formal errors). These quantities, based 
on our analysis of visual binaries, are mildly larger than the 
values derived by Heacox (1995). In the latter paper the mass 
ratio distribution is derived from a much smaller, but indepen- 
dent sample. Heacox (1995) derives the mass ratio distribution 
for spectroscopic binaries with early B-type primaries, while 
in our analysis we derive the distribution from observations of 
visual binaries with B and A-type primaries. 

10.3. Background stars 

In imaging surveys it is not always clear whether a secondary is 
a physical companion or a background star. In our models for 
the association this is obviously not a problem (as we know the 
true nature of each star in our model), but in practice it is. In 
this section we discuss briefly how background stars can afi'ect 



the interpretation of the results of a visual binary survey. We 
use the prescription derived in Kouwenhoven et al. (2007) for 
the number of background stars brighter than Kg and with an 
angular separation smaUer than p, as a function of Ks and p: 

N{Ks,p) = C-W''' -Aip), (34) 

where Aip) is the enclosed area in the field of view within a ra- 
dius p,y = 0.32 ± 0.01 mag"', and C is a constant. The Ks de- 
pendency was derived by K06 using the Besangon model of the 
Galaxy (Robin et al. 2003), and the normalization constant C is 
determined using the background star study of SHT Table 1 1 
lists for each of the three visual surveys (K05, K06, and SHT) 
the expected number of background stars per field of view. We 
list the number of background stars with Ks < 12 mag and 
with A's < 18 mag for four different pointings in the Sco OB2 
region. 

For the targets in our model, we assume that A'/3 targets 
are in each of the regions US and UCL, and that NI6 targets 
are in each of the regions LCC and the Galactic plane (GP). 
After the number of background stars A^bg is determined, each 
background star is assigned a Ks magnitude randomly drawn 
from the generating distribution corresponding to Equation 34. 
Each of the A^bg background stars is assigned a random posi- 
tion in the field of one of the A'^ target stars. Finally, the an- 
gular separation and magnitude diff'erence are calculated. For 
each background star we then decide whether it would be de- 
tected in the simulated observations, i.e., whether it satisfies the 
contrast constraint. Figure 1 1 shows the results for one of our 
models. Note how well the results in this figure resemble those 
of the K05 observations in Kouwenhoven et al. (2007) in their 
Figure 8. 

As the properties of the background star population are 
well described, and the parallax and age of each Sco OB2 
member star is known, the colour and magnitude of each sec- 
ondary can be used to determine whether it is a companion 
star or a background star. This method has been used in the 
analysis of the surveys of SHT, K05, and K06. SHT con- 
sider all new companions with Ks > 12 mag, / > 13 mag, 
ox J - Ks > 1 mag as background stars. K05 consider all 
new companions with Ks > 12 mag as background stars. The 
expected number of background stars with Ks < 12 mag in 
the K05 and SHT datasets can now be calculated. Among the 
199 targets in the K05 dataset, we expect 1 -4 background stars 
with Ks < 12 mag. In their follow-up study K06 show that the 
Ks = 12 criterion correctly classifies the companions and back- 
ground stars in about 85% of the cases. They performed multi- 
colour observations of several doubtful candidate companions, 
and identified six of these as possible background stars. Among 
the 87 confirmed Sco OB2 members in the SHT dataset we ex- 
pect ~ 1 background star with Ks < 12 mag. The expected 
number of bright background stars is small, and many back- 
ground stars have been removed by magnitude criteria and the 
follow-up study of K06. The possible presence of background 
stars among the candidate binaries thus has a negligible efi'ect 
on the conclusions of our analysis. 
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Group 


/ 


b 


Nk05 


Nk05 


Nk06 


Nk06 


NsHT 


NsHT 


C (arcsec ^) 


Kg limit 






< 12 mag 


< 18 mag 


< 12 mag 


< 18 mag 


< 12 mag 


< 18 mag 




US 


352° 


20° 


0.002 


0.19 


0.001 


0.10 


0.002 


0.15 


1.93 X 10^^ 


UCL 


328° 


13° 


0.010 


0.87 


0.006 


0.47 


0.008 


0.67 


8.71 X 10-' 


LCC 


299° 


6° 


0.016 


1.35 


0.009 


0.73 


0.013 


1.05 


13.5 X 10-9 


GP 


300° 


0° 


0.057 


4.78 


0.031 


2.59 


0.044 


3.70 


48.0 X 10-9 



Table 11. The number of background stars expected per field of view for each of the three imaging surveys discussed in this paper, according 
to our model. The field of view size is 361.2 arcsec^ for K05, 196.0 arcsec^ for K06, and 280.1 arcsec^ for SHT We list the results for 
three different pointings: to the centers of the three subgroups US, UCL, and LCC, and for the intersection of LCC with the Galactic plane. 
Columns 4, 6, and 8 list the number of background stars brighter than Kg = 12 mag per field of view. Columns 5, 7, and 9 list the number of 
background stars brighter than Ks = 18 mag per field of view. The last column lists for each group the value of the normalization constant C in 
Equation 34. 



10.4. Binary fraction versus primary mass 

Throughout this paper we have assumed that most binary pa- 
rameters are mutually independent, and independent of the bi- 
nary fraction as well. We have ignored the possible dependence 
of binary fraction, semi-major axis, eccentricity, and mass ratio 
on primary mass. However, due to the relatively small range of 
primary masses studied, these dependencies do not affect our 
results significantly. 

Several recent studies have focused on the very low-mass 
(VLM) and brown dwarf (BD) binaries in Sco OB2, such as 
the survey for spectroscopic binaries in US and p Oph by 
Kurosawa et al. (2006), and the near-infrared visual imaging 
survey among 12 VLM and BD members of Sco OB 2 by Kraus 
et al. (2005). Surveys for binarity among the lowest mass stars 
and brown dwarfs find an observed binary fraction Fm ~ 20% 
(Basri & Reiners 2006; Burgasser et aL 2007). This value is a 
lower limit to the intrinsic binary fraction among these targets 
due to the presence of undetected companions. As the details 
of the binary formation mechanism are not known, we limit 
ourselves to a discussion of two possible relations between the 
binary fraction Fm{Mi) and the primary mass Mi : a Unear and 
a logarithmic dependence. 

For simplicity, we adopt a binary fraction of 100% for 
the most massive stars in Sco OB2: Fm(20 Mq) = 100% 
and the observational constraint for the binary fraction in the 
low-mass regime: Fm(0.02 Mq) > Fm(0.02 Mq) « 20%. 
Adopting a linear relation between Fm(Mi) and Mi, we find 
thatFM(Mi) > Fm(Mi) ^ 4(Mi/M0)-H 20%. This prescription 
results in a strong dependence for high-mass stars: the binary 
fraction drops from Fm(20 M©) = 100% for the earliest B stars 
to Fm(3 Mq) > 32% for the latest B stars. If we assume, for 
example, that the intrinsic binary fraction for brown dwarfs is 
50%, the binary fraction for the latest B stars is 55%, and even 
lower for the A stars. As our observations are dominated by 
targets of spectral type A (1.5 :$ Mi $ 3 Mq), the linear pre- 
scription for Fm(Mi) can be excluded (unless the dependence 
is very weak). 

A logarithmic form for Fm(Mi) results in Fm(Mi) > 
26.7 log(Mi/Mg) + 65.3%. Then stars of spectral type B have 
a binary fraction of 80 - 100%, and those of spectral type A 
have a binary fraction of 70 - 80%. Again, if we assume an 
intrinsic binary fraction for brown dwarfs of 50%, the values 
are 85 - 100% and 80-85%, respectively. The latter values are 



in better agreement with our observations, although they are 
rather low. 

A correlation between binary fraction and spectral type 
may be present for Sco OB2, although the functional relation is 
presently unknown. Our derived value for Fm, however, is un- 
likely to be strongly affected by our choice of a binary fraction 
that is independent of spectral type, due to the relatively small 
range in mass of the binaries in our analysis. Finally, we wish 
to note that the derivation of the distributions /«(«), fp(P) and 
fqiq) are practically independent of our choice of Fm(Mi). 

11. Summary and outlook 

We have recovered the properties of the intermediate-mass 
binary population in the nearby OB association Sco OB2, 
with the aim of finding the primordial binary population 
(which is present just after star formation). We have performed 
Monte Carlo simulations, and compared for each association 
model the simulated observations with the results of surveys 
for visual, spectroscopic, and astrometric binary systems in 
Sco 0B2. The main results of our study are the following: 

- The current binary fraction among A and B stars in 
Sco OB2 is at least 70% (3<t confidence). The best agree- 
ment with the observations is obtained for models with a 
binary fraction of 100%. 

- We constrain the orbital size distribution (which can be de- 
scribed using a semi-major axis distribution or a period dis- 
tribution) using the observed angular separation distribu- 
tion and the number of visual, spectroscopic, and astromet- 
ric binaries. The observations are consistent with Opik's 
law, i.e., faia) oc a^' with ya ~ -1.0, which corresponds to 
an orbital period distribution fpiP) oc P~^. The log-normal 
period distribution found by Duquennoy & Mayor (1991), 
which corresponds to a log-normal semi-major axis distri- 
bution, is consistent with the observed angular separation 
distribution and visual binary fraction (if Fm ~ 100%), but 
significantly underpredicts the number of spectroscopic bi- 
naries. The best-fitting models have intermediate proper- 
ties, i.e., a period distribution with width crp > 2.3. 

- Random pairing (RP) and primary-constrained random 
pairing (PCRP) from the mass distribution are excluded 
with high confidence. The pairing function of intermedi- 
ate mass binaries in Sco OB2 is well-described by a mass 
ratio distribution of the form fq{q) oc q^i , with -y^ « -0.4. 
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Fig. 11. The simulated K05 data for a realization of the best-fitting 
model of Sco OB2 (cf. Kouwenhoven et al. 2007, their Figure 3.8). 
The plots show the distribution of physical companions (plusses) and 
background stars (triangles). The horizontal dashed Une indicates the 
criterion used by K05 (based on the analysis of SHT) to statisti- 
cally separate companion stars and background stars. The simulated 
sample consists of 199 targets, next to which we detect 64 compan- 
ion stars and 28 background stars. Although three companions have 
> 12 mag, and two background stars have Ks < 12 mag in 
this example, the majority of the secondaries is correctly classified 
if Ks = 12 mag is used to separate companions and background stars. 



- Sco OB2 is a young OB association with a low stellar den- 
sity. Stellar and binary evolution have practically only af- 
fected the binaries with O-type components. Dynamical 
evolution has only mildly affected the binary population. 
The current binary population of Sco OB2, as described 
above, is expected to be very similar to the primordial bi- 
nary population of Sco OB2. 

Practically all intermediate mass stars in Sco OB2 are part of a 
binary or multiple system. Although the most massive stars in 
Sco OB2 have already evolved away from the main sequence, 
observations of younger populations suggest that practially all 
of these were binary or multiple (e.g. Mason et al. 1998). This 



imphes that multiplicity is a fundamental parameter in the star 
forming process. 

In this paper we have included the results of six major 
binarity surveys among Sco OB2 members. We have not in- 
cluded the smaller surveys and individual discoveries, as each 
of these has its specific, often poorly documented selection ef- 
fects. Inclusion of these will provide a better description of 
the binary population in Sco 0B2. However, the results are 
hkely to be similar to those described in this paper, as the six 
major binary surveys used contain the large majority of the 
known binaries in Sco 0B2. Although previous observations 
(e.g. Kouwenhoven et al. 2005, 2007) have shown that triple 
and higher-order systems are present in Sco OB2, we have ne- 
glected these systems here, because of the non-trivial compar- 
ison with the observations, the more complicated selection ef- 
fects, and the very small number of known higher-order multi- 
ples. 

In our recovery of the current and primordial binary popu- 
lation in Sco OB2 we have made several assumptions and sim- 
plifications. Our main results are obtained using the visual bi- 
naries, and are shown to be consistent with the available spec- 
troscopic and astrometric binaries. The modehng of the selec- 
tion effects of the spectroscopic and astrometric binaries needs 
to be improved, so that the simulated observations can be di- 
rectly compared to the observed binary parameter distributions, 
in order to accurately derive the eccentricity distribution and 
the possible correlation between the different binary parame- 
ters. 

The assumed independence of the binary parameters, as 
well as the properties of the low-mass binary population need 
to be addressed observationally. Due to selection effects a rel- 
atively small number of binaries is known among the low- 
mass members of Sco OB2, making it difficult to derive the 
properties of these. Several studies for binarity among low- 
mass (candidate) members of Sco OB2 have been performed 
recently (e.g. Bouy et al. 2006; Kurosawa et al. 2006; Kraus 
& Hillenbrand 2007). In the near future, after a further exten- 
sion of the dataset, and a thorough membership study, it may 
be possible to characterize the binary population of Sco OB2 
over the full range of spectral types. This issue can be further 
and more accurately addressed using the results of the Gaia 
space mission (Ferryman et al. 2001; Turon et al. 2005), which 
is a project of the European Space Agency, expected to be 
launched in 2011. Gaia will survey over a billion stars in our 
Galaxy and the Local Group, and will provide an enormous 
dataset of visual, eclipsing, spectroscopic, and astrometric bi- 
naries (Soderhjelm 2005). The membership and stellar content 
of nearby OB associations can be accurately determined using 
the results of Gaia, down into the brown dwarf regime. The 
Gaia dataset will be homogeneous, and its selection effects can 
therefore be modeled in detail. The available dataset of binaries 
will be larger and more complete than any other binarity survey 
in Galactic star clusters and OB associations thus far 

The current binary population in Sco OB2 is a fossil record 
of the primordial binary population, as the young age and low 
stellar density of Sco OB2 guarantees that stellar evolution has 
affected only a handful of the most massive binaries, and sug- 
gests that dynamical evolution of the binary population has 
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been modest. The latter statement needs to be verified using nu- 
merical simulations of evolving OB associations. Whether the 
effect of dynamical evolution has been negligible over the hfe- 
time of Sco OB2 depends on its initial conditions. If Sco OB2 
was born as a low-density association, similar to its present 
state, the binary population is expected to have changed only 
modestly due to dynamical evolution. On the other hand, if the 
association has expanded significantly over the last 5-20 Myr, 
dynamical evolution may have been more prominent. 
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Appendix A: Datasets used 



Table A.l: Properties of the K05 dataset that we use in our analysis. Note that not all candidate companions are listed here, as 
in our analysis we only consider single stars and binary systems. The members HIP68532 and HIP69113 (marked with a star) 
both have two companions with a similar separation, position angle, and brightness. In our analysis we consider these "double 
companions" as a single companion. 
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HIP75915 


6.44 


8.15 


1.71 


5.60 


229 


2.89 


1.22 


0.42 


HIP76001 


7.60 


7.80 


0.20 


0.25 


3 


1.54 


1.36 


0.88 


HIP76071 


7.06 


10.87 


3.81 


0.69 


41 


2.70 


0.23 


0.09 


HIP77315 


7.24 


7.92 


0.68 


0.68 


67 


2.08 


1.56 


0.75 


HIP77911 


6.68 


11.84 


5.16 


7.96 


279 


2.80 


0.09 


0.03 


HIP77939 


6.56 


8.09 


1.53 


0.52 


119 


3.85 


1.82 


0.47 


HIP78756 


7.16 


9.52 


2.36 


8.63 


216 


2.30 


0.92 


0.40 


HIP78809 


7.51 


10.26 


2.75 


1.18 


26 


2.03 


0.30 


0.15 


HIP78847 


7.32 


11.30 


3.98 


8.95 


164 


2.20 


0.160 


0.07 


HIP78853 


7.50 


8.45 


0.95 


1.99 


270 


1.82 


1.14 


0.63 


HIP78956 


7.57 


9.04 


1.47 


1.02 


49 


2.40 


1.16 


0.48 


HIP79124 


7.13 


10.38 


3.25 


1.02 


96 


2.48 


0.33 


0.13 


HIP79156 


7.61 


10.77 


3.16 


0.89 


59 


2.09 


0.27 


0.13 


HIP79250 


7.49 


10.71 


3.22 


0.62 


181 


1.42 


0.140 


0.10 


HIP79530 


6.60 


8.34 


1.74 


1.69 


220 


3.73 


1.58 


0.42 


HIP79631 


7.17 


7.61 


0.44 


2.94 


128 


1.90 


1.58 


0.83 


Continued on next page 



Kouwenhoven et al.: Recovering the binary population for intermediate mass stars in Sco 0B2, Online Material p 3 



Table A.l - continued from previous page 



fflp 


Ks,i 


Ks,2 


AKs 


P 


<f 


Ml 


M2 


q 




mag 


mag 


mag 


arcsec 


deg 


Mo 


Mo 




fflP79739 


7.08 


11.23 


4.15 


0.96 


118 


2.32 


0.16 


0.07 


HIP79771 


7.10 


10.89 


3.79 


3.67 


313 


2.14 


0.19 


0.09 


HIP80238 


7.34 


7.49 


0.15 


1.03 


318 


1.94 


1.67 


0.86 


HIP80324 


7.33 


7.52 


0.19 


6.23 


152 


1.70 


1.54 


0.91 


HIP80371 


6.40 


8.92 


2.52 


2.73 


141 


3.43 


0.94 


0.27 


HIP80425 


7.40 


8.63 


1.23 


0.60 


156 


2.08 


1.16 


0.56 


fflP80461 


5.92 


7.09 


1.17 


0.27 


286 


5.29 


2.97 


0.56 


HIP80799 


7.45 


9.80 


2.35 


2.94 


205 


1.86 


0.34 


0.18 


HIP80896 


7.44 


10.33 


2.89 


2.28 


177 


1.81 


0.24 


0.13 


HIP81624 


5.80 


7.95 


2.15 


1.13 


224 


6.53 


2.30 


0.35 


HIPS 1972 


5.87 


11.77 


5.90 


5.04 


213 


4.92 


0.35 


0.07 


HIP83542 


5.38 


9.90 


4.52 


8.86 


196 


1.10 


0.91 


0.83 


HIP83693 


5.69 


9.26 


3.57 


5.82 


78 


4.95 


1.06 


0.21 



Table A.2: Properties of the K06 dataset that we use in our analysis. Note that not all candidate companions are Usted here, as 
in our analysis we only consider single stars and binary systems. The members HIP68532 and HIP69113 (marked with a star) 
both have two companions with a similar separation, position angle, and brightness. In our analysis we consider these "double 
companions" as a single companion. 



HIP 


KsA 


Ks.2 


AKs 


P 


V 


M, 


M2 


q 




mag 


mag 


mag 


arcsec 


deg 


Me 


Ms 




H1P59502 


6.87 


11.64 


4.77 


2.935 


26 


1.80 


0.14 


0.08 


HIP60851 


6.06 


13.69 


7.63 


8.159 


231 


2.63 


0.04 


0.02 


HIP61265 


7.46 


11.38 


3.93 


2.505 


67 


1.82 


0.27 


0.15 


HIP62026 


6.31 


7.86 


1.55 


0.232 


6 


2.45 


1.19 


0.49 


fflP63204 


6.78 


8.40 


1.62 


0.153 


237 


2.05 


1.06 


0.52 


HIP67260 


6.98 


8.36 


1.38 


0.423 


229 


2.00 


1.10 


0.55 


HIP67919 


6.59 


9.10 


2.51 


0.685 


297 


1.97 


0.75 


0.38 


HIP68532* 


7.02 


9.54 


2.53 


3.052 


288 


1.95 


1.12 


0.57 


HIP69113* 


6.37 


10.29 


3.92 


5.344 


65 


3.87 


1.49 


0.39 


HIP73937 


6.23 


8.37 


2.14 


0.242 


191 


2.94 


1.11 


0.38 


HIP78968 


7.42 


14.26 


6.84 


2.776 


322 


2.33 


0.02 


0.01 


HIP79739 


7.08 


11.23 


4.15 


0.959 


118 


2.32 


0.16 


0.07 


HIP79771 


7.10 


11.42 


4.33 


0.435 


129 


2.14 


0.19 


0.09 


HIP80799 


7.45 


9.80 


2.35 


2.940 


205 


1.86 


0.34 


0.18 


HIP80896 


7.44 


10.33 


2.89 


2.278 


177 


1.81 


0.24 


0.13 


HIP81949 


7.33 


15.52 


8.19 


5.269 


341 


2.26 


0.02 


0.01 


HIP81972 


5.87 


11.77 


5.90 


5.040 


213 


4.92 


0.35 


0.07 


HIP83542 


5.38 


9.90 


4.52 


8.864 


196 


1.10 


0.91 


0.83 
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Table A.3. Properties of the SHT dataset that we use in our analysis. Note that not all candidate companions are listed here, as in our analysis 

we only consider single stars and binary systems. The six members at the bottom of the list were not explicitly observed by SHT. Due to the 
presence of (known) close companions these were not suitable for wavefront sensing. We have included these targets for our analysis to avoid 
a bias towards low binarity. 



HIP 


Ks,i 


Ks,2 


AKs 


P 


<P 


Ml 


M2 






mag 


mag 


mag 


arc sec 


deg 


Mo 


Mo 




H1P55425 


4.66 


5.86 


1.20 


0.354 


144 


4.65 


2.70 


0.58 


HIP56561 


3.17 


6.81 


3.64 


0.734 


135 


8.30 


2.22 


0.27 


H1P58884 


5.67 


7.00 


1.33 


0.698 


158 


3.17 


1.75 


0.55 


H1P61585 


3.41 


10.94 


7.53 


4.853 


198 


6.30 


0.19 


0.03 


H1P63945 


5.80 


9.16 


3.36 


1.551 


268 


3.60 


0.135 


0.04 


HIP65271 


5.12 


7.03 


1.91 


0.164 


135 


4.25 


1.80 


0.42 


HIP67472 


3.97 


10.06 


6.09 


4.637 


304 


7.95 


0.75 


0.09 


H1P72683 


5.27 


6.84 


1.57 


0.099 


86 


4.52 


2.17 


0.48 


H1P72800 


5.54 


9.43 


3.89 


1.046 


161 


3.75 


0.76 


0.20 


rlir / jjj^- 


A no 




1 '^7 


U. IZo 




/ .OJ 




U.Oo 


HIP75264 


4.28 


5.55 


1.27 


0.279 


149 


7.25 


4.84 


0.67 


HIP76945 


5.79 


9.47 


3.68 


0.507 


133 


3.50 


0.80 


0.23 


H1P77939 


6.13 


7.78 


1.65 


0.524 


120 


4.95 


2.16 


0.44 


H1P78820 


3.86 


6.80 


2.94 


0.292 


171 


11.20 


2.98 


0.27 


HIP79374 


4.20 


5.14 


0.94 


1.334 


2 


8.32 


5.47 


0.66 


HIP79530 


6.31 


8.07 


1.76 


1.693 


220 


3.23 


1.35 


0.42 


HIP80112 


2.61 


4.77 


2.16 


0.469 


244 


19.96 


10.40 


0.52 


H1P57851 








1.549 


158 


4.15 


1.83 


0.44 


HIP62322 








1.206 


35 


7.35 


6.40 


0.87 


HIP64425 








0.185 


7 


4.07 


3.10 


0.76 


H1P74117 








0.193 


210 


6.49 


4.94 


0.76 


H1P76371 








2.150 


8 


5.75 


2.79 


0.49 


HIP77840 








2.162 


270 


6.05 


2.39 


0.40 



Table A.4. The LEV dataset used in this sample, consisting of 16 binaries with orbital elements (SBl or SB2), and 23 radial velocity variables 
(RVV), for which no orbital elements are available. Left: the 16 spectroscopic binaries with orbital elements among the confirmed members of 
Sco OB2, in the LEV dataset. LEV observed 53 confirmed members of Sco 0B2, of which 8 SB Is, 8 SB2s, 23 RVVs, and 14 targets with a 
constant radial velocity. Right: the 23 RVVs. Note that several of the RVVs may not be spectroscopic binaries, as radial velocity variation may 
also be caused by line profile variabiUty. 



HIP 


P 


e 


a> 


Group 


Radial velocity variables 




days 




deg 








HIP67464 


2.6253 


0.13 


222 


UCL 


HIP67472 


HIP68245 


HIP75647 


3.8275 


0.25 


22 


UCL 


HIP68862 


HIP70300 


HIP76297 


2.8081 


0.10 


97 


ULC 


HIP73334 


HIP74100 


HIP76503 


5.2766 


0.33 


86 


US 


HIP75141 


HIP75304 


HIP76600 


3.2907 


0.28 


114 


UCL 


H1P76633 


H1P77635 


HIP76945 


12.26 


0.19 


83 


UCL 


HIP77900 


HIP77939 


HIP77858 


1.9235 


0.36 


309 


US 


HIP78246 


HIP78384 


HIP77911 


1.264 


0.61 


330 


us 


HIP78530 


HIP78549 


HIP78104 


4.0031 


0.27 


231 


us 


HIP78655 


HIP79031 


HIP78168 


10.0535 


0.58 


340 


us 


H1P79374 


H1P79739 


HIP78265 


1.5701 


0.15 


25 


us 


HIP80024 


HIP81266 


HIP78820 


6.8281 


0.28 


38 


us 


HIP82545 




HIP79404 


5.7805 


0.19 


115 


us 






H1P79374 


5.5521 


0.11 


267 


us 






HIP80112 


34.23 


0.36 


308 


us 






HIP80569 


138.8 


0.44 


325 


us 







